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Abstract 

For classical Hamiltonian A^-body systems with mildly regular pair in- 
teraction potential (in particular, Hf^^ integrability is required) it is 
shown that when ^ oo in a fixed bounded domain A C M^, with 
energy £ scaling as £ oc A^^, then Boltzmann's ergodic ensemble en- 
tropy Sa{N, £) has the asymptotic expansion S\{N, N'^e) = — A^ln A^-|- 
s^{e)N + o{N); here, the A^lnA^ term is combinatorial in origin and 
independent of the rescaled Hamiltonian, while s^C^) system- 
specific Boltzmann entropy per particle, i.e. — s^(e) is the minimum 
of Boltzmann's H function for a perfect gas of energy e subjected to a 
combination of externally and self-generated fields. It is also shown that 
any limit point of the n-point marginal ensemble measures is a linear 
convex superposition of n-fold products of the //-function-minimizing 
one-point functions. The proofs are direct, in the sense that (a) the 
map £ I-)- 5(£) is studied rather than its inverse S i-7> -E(S); (b) no reg- 
ularization of the microcanonical measure 5{£. — H) is invoked, and (c) 
no detour via the canonical ensemble. The proofs hold irrespective of 
whether microcanonical and canonical ensembles are equivalent or not. 
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1 Introduction 



The rigorous foundations of equilibrium statistical mechanics have largely been 
laid long ago |Rue69t IPenTOt ILen73t IML79] . but the most basic problem in 
classical statistical mechanics, namely the rigorous asymptotic evaluation of 
Gibbs' microcanonical ensemble |Gib02] in the limit of a large number 
of particles, has only been treated in an approximate way. The standard 
way of dealing with the microcanonical ensemble (a.k.a. Boltzmann's ergodic 
ensemble |Bol96] ) in a rigorous manner |Rue69t ILan73t IML79] has been to 
replace its singular ensemble measure by a regularized measure (usually also 
referred to as microcanonical, although quasi-microcanonical would seem a 
better name). In these approaches one cannot take the limit of vanishing 
regularization; yet, since one can approximate the singular measure as closely 
as one pleases, "this is not completely unsatisfactory from a conceptual point 
of view" ( |Lan73] . p. 4). All the same, Lanford's wording makes it plain that it 
is desirable to find a way to remove the regularization or to avoid it altogether. 

Recently [Kie09aj the author noticed that after only minor modifications, 
Ruelle's method |Rue69] to establish the thermodynamic limit for Boltzmann's 
ergodic ensemble entropy, taken per volume (or per particle), works without 
the need for any regularization of the ensemble measure; a follow-up work on 
the thermodynamic limit of the correlation functions is planned. Taking "t/ie 
thermodynamic limit" |Rue69j means that the domain A grows "evenly" with 
A^ and such that N/Vo\ (A) — >■ p with p a fixed number density, and the energy 
£ scales such that £/Vol (A) — > e (or 8,/N ^ e, abusing notation), with e G M 
a fixed energy density (or energy per particle) — this limit covers systems of 
interest in condensed matter physics or chemical physics, such as those with 
hard core or Lennard- Jones interactions. 

In the present paper we will be concerned with another limit A^ — > oo, 
where A is fixed and £ scales such that £/A^^ — )■ e. This limit covers sys- 
tems of interest in plasma and astrophysics, such as those with Coulomb or 
(mollified) Newton interactions. It is variably known0 as a "thermodynamic 
mean-field limit," a "self-averaging limit," or "Vlasov limit." We will study 
the Boltzmann ensemble entropy and the correlation functions. 

The remainder of this paper is structured as follows. In section 2 we collect 
the defining formulas of the ergodic / microcanonical ensemble for finite A^ and 
explain which probabilistic quantities are of physical interest. In section 3 we 
give a heuristic motivation for the Vlasov limit. In section 4 we state our main 
theorems, ordered by increasing depth. Their proofs are given in sections 5.1 
to 5.3. Section 6 lists some spin-offs of our results, and section 7 closes our 
paper with an outlook on some open problems. 

^Thc first two names refer to a Weiss- type "mean-field approximation" becoming exact 
in the fimit, but we will not invoke any such approximation and speak of the Vlasov limit. 
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2 A brief review of the ergodic ensemble 

For a Newtonian A^-body systenH in a domain A C M'^ with Hamiltonianlfl 

l<j<Af l<i<j<N l<i<N 

the ergodic / microcanonical ensemble is a family {x[,^''|A; G N} of i.i.d. copies 
of a random vector X^^) = (Pi, Qi; . . . ; P^v, Qtv) e (M^ x A)^ distributed 
according to the stationarjfl single-system a-priori probability measure 

A 

where := (Pi, Qi; Pat, ^at) G (M^xA)^ and d^^X is 6X-dimensional 

Lebesgue measure, and where 

fi;(^)(e) = hj^{^- Hi''\x(^^)) d'^X (3) 

is known as the structure function^ here, the ' means derivative w.r.t. £ of 

= W.J X{Hf)<£}d'^^' (4) 

where X^h^n)^^^ is the characteristic function of the set {H^\x'-^'>) < £} G 
(M^xA)^, over which the integrals extend. Thus, if ^ denotes the Borel sets 
of (M^xA)^ C M^^, then ((M^x A)^, ^, /ij.^^) is the single-system probability 
space; so if i? G is a Borel set, then the probability of X^^^^ being in B i^ 

Prob(X(^) G 5) =/if^(B). (5) 

Clearly, Prob (X^^) G fi) = unless B f] {i/f ^ = £} ^ 0; put differently, 

Fioh{H'^\x.[^^) = £) = 1 V A; G N. Moreover, Prob (X(^) G d^^X) = 

lj,[^\d^^X) is the a-priori probability for X(^) to be in d^^X about X(^). 

The ergodic ensemble is probabilistically meaningful for all X G N, yet its 
thermodynamic significance emerges only in the large X regime (Avogadro's 
X 10^^) when it makes sense to speak of a solid, a liquid, a plasma (etc.) on 
macroscopic scales of space and time. Since the typical physical characteristics 



^ All particles belong to a single specie. We use units of mc^ for energy, mc for momentum, 
and h/mc for length, where m is particle mass, c the speed of light, h Planck's constant. 

■^It is understood that Wa is symmetric, i.e. WA{q,q) — WA{q,q), and not reducible to 
a sum of one-body terms; other details of W and V will be specified in the next section. 

"^Stationarity is defined w.r.t. the flow generated by the Hamiltonian H^\pi, . . . , g^). 

^The A^! term was supplied by Gibbs to resolve Gibbs' paradox. It cancels out in ([2]). 

^It is tacitly understood that whenever one encounters a physically interesting subset L 
of a Borel null-set which is not itself Borel measurable, then we use the Lebesgue cr-algebra. 
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of solids and liquids (etc.) are not revealed by "picturing" such systems as 
individual points in R^^, one associates each microstate X^^^ with a unique 
family of empirical n-point "densities" on (M^xA)", n = 1,2,. ..,N. The 
normalized one-point "density" with N atoms (empirical measure) is given by 

A^;!,, (p, q) = ^ 6ip- - Qi) (6) 

l<i<Af 

and the normalized two-point density with N atoms ([/-statistic of order 2) by 

A^Mq-^p'^q') = iv(F3T) EE'^(p-P^)'^(9-g.)5(p'-P,)'5(g'-9,); (7) 

l<i^j<N 

similarly the empirical n-point densities with n = 3,...,N are defined. The 
map X^^) — > {A^^lM)}n=i is bijective if we insist that the particular labeling 
given to us algebraically with r.h.s.dS]) or r.h.s.(I7]) etc. has an intrinsic mean- 
ing; however, considered purely measure theoretically as "density" on M®" each 
A^(^) is invariant under the permutation group applied to the particular label- 
ing, and since there are N\ distinct X^^^s obtained by permuting the particle 
labels, the map X^'^^ — >■ {A^(jv)}^=i is many-to-one in this sense. Understood 
in this measure theoretic way the empirical ra-point densities do not depend on 
the unphysical (though mathematically convenient) labeling of the particlesjll 
and so are physically more natural than points in R^^; when n is small, say 
n = 1 or 2, then A^jj^, is also "physically more manifest" than a point in 
Hence, the probabilities of interest to physicists will be of the form 

Prob (a JJ^, e B 

for physically significant measurable sets B in *P'^((R^xA)"), the permutation- 
symmetric probability measures on (R^xA)". Among the physically significant 
sets are balls (w.r.t. a suitable topology, still to be chosen) centered at a rep- 
resentative n-point density function for a solid, liquid, ... , or complements of 
such balls. As for the topology, the fine {TV) topology for *^'*((R^ x A)") is 
not suitable as it is equivalent to discriminating between different ^^\n) w.r.t. 
the Borel sigma algebra of R^^/Sat (see footnote 6). Practically accessi 
are only some considerably less finely resolved events, such as the empirical 
n-point densities ^^(n) distinguished w.r.t. the weak topology, quantified by 



^So, physically we can identify these A^! distinct X^^^s with a single N-point configuration 
m R'^ X A, which is a point X(^) G ^^^xK^/Sn. The subscript ^ means that coincidence 

points are removed, and Sn is the symmetric group of order N. We should also write A^-'^^ , 
with the understanding that as measure A*^\n) given by A^^, for any of the A^! points 
X(^) in the pre-image in of X'^^\ The map X^^) ^ W~\^-,}n^i is bijective. 

^Even if the balls in TV topology were practically accessible, for iV 3> 1 the amount of 
information would be sheer overwhelming and not very illuminating. 
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a convenient Kantorovich-Rubinstein metric c^kr on ^*((M^xA)"). Two very 
different points in R^/Sn, say X(^) and Y^^\ can map into two densities 
A^(^) and ^yi^) which in weak topology on *P*((]R'^xA)"') are virtually indis- 
tinguishable; here X^^^ and Y^^^ are any two representative points out of the 
N\ points each, which constitute the pre-image in of X^'^\ respectively 
Y^^\ So the probabilities of interest to physicists are typically of the form 

FToh{dKK{A^^l,j!:^^)>6y (9) 

where 9^^^; g^")) G (q3n(£g)((M3xA)") is an equilibrium density 

function, defined — in the simplest of all cases — implicitly as the unique 
function for which (after rescaling of variables and parameters, if necessary) 

Prob (^KR (aJ),) , /i?) > 5) "-^ V5 > 0. (10) 

In these simplest of all cases, ( JTOj) also explains what is meant by a "representa- 
tive n-point density function;" and whether /iq^ represents a solid, liquid, gas, 
etc., depends on the specific configurational correlations exhibited by /iq\ 
In more complicated (and more interesting) situations, several "competing" 
equilibrium functions /iq^ may exist, and ( ITOl) has to be modified accordingly. 

The "simplest case" scenario just described was discovered by Boltzmann 
(p. 442 of |Bol96] ). based on his explicit evaluation of (|2]) for the perfect gas. 
He realized that when H^'' is the perfect gas Hamiltonian and ^ 1, then 
basically every point of {H^^ = £} (identified with an n-pt. density through 
the map X^^'^ — )■ A^^^^j) lies in the vicinity (w.r.t. weak topology) of one 
and the same equilibrium density function fiq^ at that energy £, and given n. 
When sports non-trivial pair interactions, Boltzmann's description needs 
to be modified slightly to account for the phenomenon of phase transitions. 

While there can hardly be a doubt that Boltzmann's insight into ([2]) is 
correct, the rigorous results which support his assessment have been obtained 
not for ([2]) but for some regularized approximation of this singular measure 
|Rue69t ILen73t IML79] . In this paper we will finally vindicate Boltzmann's 
ideas in the Vlasov regime of the relevant class of Hamiltonians ([1]) . 

3 Heuristic considerations on the Vlasov limit 

For the ergodic ensemble to exhibit a Vlasov regime the Hamiltonian ([1]) 
needs to satisfy additional conditions. In particular, a necessary condition 
on the symmetric and irreducible pair potential W\ is local integrability, i.e. 
Wa{Q} ■) ^ ^^{Br{q) n A) Vq e A. We remark that for the existence of a 
dynamical Vlasov regime the local integrability of the forces derived from Wa 
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is mandatory, viz. VqWaIq, ■ ) G £,^{Br{q) fl A) Vq G A. Coulomb's electri- 
cal and Newton's gravitational interactions belong in either class. Physically 
meaningful external potentials V^^"* are continuous for q G A; it has minor 
technical advantages to assume that V^^'' is actually continuous also at the 
boundary, i.e. \imqi^qVl^\q') = vl^\q) for all q E dA and q' G A. For 
convenience we assume that inf Hj^\pi, . . . , q^) = min Hj^\pi^ . . , q^y) = 
^g{N) > — oo, and call £g(iV) the A^-body ground state energyp Newton's 
gravitational interactions need to be regularized to achieve 8,g{N) > — oo. 

In the introduction we have already mentioned that the Vlasov limit scaling 
for such interactions is £ x N'^e for N ^ 1. We now explain why. Integrating 
(E]) over p-space gives a normalized one-point "density" (empirical measure) 
on A with N atoms, which by abuse of notation we denote as follows. 



Whenever Boltzmann's simplest scenario holds, then there is an equilibrium 
density ps^at G (^n(£°)(A), depending on A^(^ 1) and £, such that A^|jv) (q) ~ 
P£.,N{q) for overwhelmingly most X^^^ distributed by (E]), where "~" means 
the two "densities" do not differ by much in a conventional Kantorovich- 
Rubinstein metric c^kr- This suggests that when A C is fixed and N ^ oo 
together with £ — > oo such that 8./N°' — )■ e for a yet-to-be determined a, 



then p£,;v p. G (^ n £0)(A) and A^iU{q)d^p Pe(q), weakly. Im- 



plementing this law-of-large-numbers type scenario inevitably leads to a = 2, 
as is most easily seen if we assume for a moment that G (t^(AxA). Then 
q (-)■ W^iq, q) is a bounded continuous function in A and we can write 



^Presumably boundedness below is not technically necessary. We expect that pair inter- 
actions which diverge logarithmically to — cx) can be accommodated but require additional 
weak compactness estimates, e.g. in some 2P space; cf. |KiLe97) . 






(12) 




q)^l\N) (P, g)d V'gA^wiv) (P, q)d3pd=^g, 
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and when J^3 A^|iv) (P, q)d^P ~ Pe{q), we find 

^(iv)(XW) ^nJJ i|ppA«,,(p,9)dV='g 

+iV 1(^/^(9) - iW^A(q,g))p,(g)d3g (13) 



+1 



The last term clearly scales oc A^^ because Wa and are independent of A^. 
In a sense this already establishes the £ oc A^^ scaling. However, we have yet 
to consider the terms on the first two lines on the r.h.s. of f|T3l) . It would 
seem that these scale oc A^ and so, for large A^, would become insignificant 
as compared to the one in the last line, but only the A^ / ^W\{q, q)p!r{q)d^q 
contribution will surely become insig nifican10 for large A^, for the same reasons 
for why the last one scales oc A^^ {Wa and do not depend on A^). As for the 
external potential vl^\q), the superscript indicates that we may want to 
adjust it to the number of particles in the system on which it acts in order to 
retain a noticeable effect when A^ becomes large. So in particular we can set 
Vi'^\q) = NVAiq) [or = (A^ - l)^A(g)], with VA{q) independent of A^, and 
find A^ Jvl^\q)p,{q)d^q = J VA{q)pe{q)d^q [+0{N)], scaling oc A^^ [j^ 
leading order], hence remaining significant in f[T^ as A^ becomes large. And as 
to the kinetic energy term, it is important to realize that J^^ A^liv) (P; q)d'^p ~ 
Peiq) e (^ n C:0)(M3) does not imply that A^^l^^{p,q) ^ f,{p,q) G (^3 n 

(r^)(]R3xA). For instance, we can have that A^^/^ A^j^, (A^^/ V, q) ~ fe{p, O) e 
(*pn £°)(M^xA) so that a significant fraction of the energy will be distributed 
over the kinetic degrees of freedom]^ and then, up to terms of 0{N), we find 

^(A^)(XW) ^ n2 (^jj {l\p\^ + VA{q)) fe{p,q)dM''q (14) 

WA{q, q)fs{P, q)fe{P, q)d^pd^qd^pd^q 




This scaling scenario can be verified explicitly for the perfect gas (Wa = 0) by 
inspecting Boltzmann's calculations, and it is reasonable to expect that it will 
continue to hold for a physically interesting class of Wa ^ 0. 

To summarize, the Vlasov limit for the Hamiltonian ([T]) with V^^^ = NV 

means that N^/^A'^\^^{N^/^p,q) feip,q) weakly in ^^(M^ x A), with 



^^Incidentally, this indicates that the Vlasov limit does not require the continuity of W\, 
the only purpose of which was to furnish identity (fT2|) which involves WA{q, q)- 

-'^-'^Unless £ is the ground state energy for which all particle momenta vanish, indeed. 
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fe{p, q) G n C:0)(M3x A), and Ar-2i7W(xW) £(/,) = e > e^, where 

+ nil ^WA{q, q)f{p, q)f{p, ~q)d'pd'qd'pd'q (15) 
is the "energy of /," and where Eg = inf /6fp(R3xA) £ (/) is given by 

Sg = Jnf^^ ( I V^{q)p{q)d\ + ^ \W^{q, ~q) p{q) p{q)d\d^) . (16) 

4 The Vlasov limit for Boltzmann's Ergode 

We now state our main results about the Vlasov scaling limit for Boltzmann's 
ergodic ensemble of A^-body systems in a format which will be recognized as 
the familiar folklore by anyone with a joint expertise in Vlasov theory and 
statistical mechanics. We will also utilize some less familiar notions. 

In the following, A C M'^ is a bounded, connected domain (open) which 
does not depend on A^. The upshot of the previous section is that if we want 
the external potential to remain significant when A^ gets large, then our A^- 
body dynamics in A will be governed by Hamiltonians ([1]) of the special type 

l<j<Af l<i<j<N 

with the single particle potential Va and the pair interaction W\ independent 
of N . We choose — 1 rather than A^ as scaling for V^^"* because then we can 
absorb Va and W\ together in a new A^-independent effective pair interaction 
Uhiq, q) '■= W\{q, q) + VA(q) + V\{q). This doesn't affect any of our results 
(as we will prove), but the Hamiltonian (IT7|) can be recast shorter as 

i/f)(Pi,...,g^)= + (18) 

l<i<N l<i<3<N 



Since heuristically we expect for a Hamiltonian system with Hamiltonian (|TI 
under Vlasov scaling that A^^''^ A^^^, (A^^/^p, g) ^-^^ fe{p,q) weakly, with 
fe{P: q) G (^n(t°)(]R^xA), we also expect that the rescaled particle momentum 
random vectors A^~^/^Pj converge in distribution, implying that | |Pi|^ ~ 
A^^^kin for /i£^''-most X*^^\ where ^kin is the kinetic energy contribution to e. 
We find it more convenient to work with random variables which themselves 
converge in distribution and so re-scale the momentum variables as p^. = 
N^/'^Pj^ in (ITS!) : or in more economical notation: we replace — )■ N^^'^p^ in 
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(fT8|) . With this minor additional abuse of notation our Hamiltonian finally 
reads 



(19) 



l<i<N l<i<j<N 



Our main results will be proved under the following hypotheses on U\{q, q): 

(HI) Symmetry: U^iq, q) = UA{q, q) 

(if 2) Lower Semi- Continuity. UA{q,q) is l.s.c. on Ax A 



^ ^ ^' JJ '^{UA{q,q)-mmUA<e\ ^ ^ 

(HA) Local Square Integrability: [^(q, ■) E S."^ {Br{q) HA) V q G A 
(ii5) Confinement: Uxiq, q) = +00 whenever q ^ A or q ^ A 

Hypothesis (HI) is a consequence for W\ of Newton's "actio equals re-actio," 
plus the symmetrized added contribution of Va, both of which need no further 
commentary. Hypothesis {H2) is satisfied by many important pair interac- 
tions invoked in physics, though not by all. For instance, the Coulomb pair 
potential U^°^\q,q) = l/\q — q\ for q ^ q satisfies {H2) after also setting 
Ul°^\q,q) = u for any particular m G M. On the other hand, the New- 
ton pair potential Uj^'^'^^{q, q) = —Ul"'^\q,q) does not satisfy {H2) for any 
choice of u; however, the regularized Newton pair potential U^^'^g{q,q) = 
^(Xfl,, * f^A°"' * Xi3^)(9) q) (where f * g denotes the conventional convolution 
product of / and g) does satisfy {H2). By {H2), there exists an A^-dependent 
ground state energy £,g{N), i.e. > 8,g{N) > —00, but the ground state 

configuration can have some unwanted features]^ Hypothesis {H3) eliminates 
the possibility of energetically isolated ground states, thus guaranteeing the 
existence of a fat set of minimizing sequences of configurations. Hypothesis 
[HA) is a little stronger than necessary, but it allows us to make convenient 
use of Chebychev's inequality to prove a law of large numbers for the pair- 
specific interaction energy; the important Coulomb potential satisfies {H4). 
Note that (-^4) implies local integrability of Ua, which is needed in var- 
ious integrals featuring in the Vlasov limit. Note also that by {H2)&i{H4) 
there exists an A^- independent G M defined by (fT6!) . In Appendix A we 
show that {H1)^{H2) guarantee that the pair-specific ground state energy 
Eg{N)/[N{N — 1)] = eg{N) is monotonic increasing with A^, and using also 
{H3) and (HA) we show that £g{N) Eg as N 00. In Appendix A we 
also show that if Ua > 0, then also tg{N)/N'^ = eg{N) Eg as N ^ 00. 

^^For instance, in our example of the amended Coulomb pair potential one can choose 
u = 0, but then Thomson's problem on §^ C [Tho04| yields as ground state configuration 
always the spurious one (up to SO{3) action) for which all particle positions coincide. To 
avoid these spurious ground state configurations it is advisable to choose u > huge. 
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Hypothesis {H5) is really inherited from the dynamical theory of particles 
in A C M^, where one sets V^^"* = +00 for qr ^ A to dynamically model con- 
finement in a container; {H5) has a minor notational advantage by allowing 
us to treat physical space integrals like momentum space integrals as over all 
M^, the spatial cutoff to A automatically being provided by the potential Va 
through [/a- Usually, {H5) is not listed explicitly as a hypothesis on the in- 
teractions even when spatial integrations are explicitly restricted to A. This 
concludes our commentary on the list of hypotheses (HI) — {H5). 

All our results (except Proposition[7]in Appendix A) will be formulated and 

proved under the convenient assumption that Ua > 0, so that Eg > 0. Since 

2 

Ua has a minimum in A , by {H2), and since the physics of our dynamical 
system does not change if we simply add a constant to Ua, we may assume that 
Ua> without loss of generality. We emphasize that this choice is merely for 
convenience, given {H2), and so is not listed as another hypothesis. 

The simplest objects of interest are the thermodynamic functions. In the 
1960s and hence, techniques based on monotonicity, convexity and super- 
additivity estimates have been developed to prove their existence and regular- 
ity in the limit N 00 which avoids having to control the more sophisticated 
objects of interest, which are the correlation functions. For the traditional 
thermodynamic limit scaling, see Ruelle's book |Rue69 ] and |Kie09a] for a re- 
cent extension of Ruelle's arguments to Boltzmann's Ergode proper. For the 
Vlasov scaling of the canonical ensemble, see |Kie93] . To extend these argu- 
ments to Boltzmann's Ergode proper with Vlasov scaling, our first goal is to 
show that the logarithm of the structure function (|3l) for the Hamiltonian (|T9|) . 
which yields Boltzmann's ergodic ensemble entropjtj (cf. eq.(305) in |Gib02] ). 

5^(^)(£)=lnn' (£), (20) 

admits the correct type of asymptotic expansion for N 00 with £ = N'^e, 
and has the correct quahtative e dependence. The usual strategy can be put 
to work if we assume just a little more than {H1)-{H5). In this vein we state: 

Theorem 1. Let H^^^ be given in / fTPj) . with Ua satisfying conditions {HI) 
and {H5), but with {H2), {H3), {HA) replaced by the single stronger condition: 

{H6) Continuity: UA{q,q) is continuous on Ax A. (21) 

Let e > Eg, with Eg > defined as before. Then the ergodic ensemble entropy 
has the following asymptotic expansion for N ^ 1, 

S^w {Nh) = -N\nN + Ns^ {e) + o{N) , (22) 

where ^^(e) is the system- specific Boltzmann entropy per particle. The function 
e (-)■ Sf^{£) is continuous and strictly increasing for e > Eg. 

^■^Entropy is measured in units of /cb, where fee is Boltzmann's constant. 
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We remark that the leading term of r.h.s. (l22| ) is purely combinatorial in 
origin and independent of the Hamiltonian — it is solely due to the A^! 
in (|3]). System-specific information begins to show in the next to leading term, 
which is 0{N). The o{N) term in ( !22l) is presumably O(lnA^). 

We will also prove two upgrades of Theorem 1 (Theorems 1^ and 1^^) 
which involve the decomposition of the system-specific Boltzmann entropy 
per particle s^{e) into a "kinetic" and an "interaction" contribution. The 
discussion of this more technical material is postponed until section 5.1. 

While they do yield valuable qualitative information about the thermo- 
dynamic functions for the systems under study, in this case s^{e), existence 
theorems such as Theorem [1] and their "proofs by sub-additivity" have the 
disadvantage that they do not characterize the limit objects in a way which 
would allow their systematic evaluation for physically interesting irreducible 
pair potentials W\ and external one-body potentials Va- It is this type of 
characterization that we are after, and in section 5.2 we prove that s^{e) sat- 
isfies the familiar maximum entropy variational principle for the entropy per 
particle of a perfect gas in a combination of self- and externally generated 
fields. More precisely, we prove the following strengthening of Theorem [1] 

Theorem 2. Let h';^^ be given in Iil9\). with Ua > satisfying (H1)-(H5) . 
Let e > Eg. Then the Boltzmann entropy [2D^) has the asymptotic expansion 

Smn) {N^e) = -N\nN + Ns^ {e) + o{N) (23) 

A 

for N ^ 1, and the system- specific Boltzmann entropy per particle is given by 
s,{s) = -Mfe), (24) 
where ^b(/) is "Boltzmann's H function" of f, which read^ 

Mf) = II fiP, q) HfiP, q)/e)d'pd% (25) 

and where f^ is any minimizer of this H functional over the set of trial densities 
21^ = {/ G (<pn£i n£Mn£i)(R3xA) ■.■E{f)= e}, where £(/) now reads 

£ (/) = ll\\P?f{P. q)dVq+llll |t/A(q, ~q)f{p, q)f{p, q)d'pd'qd'pd% (26) 
Any minimizer f^ of ^sif) over the set 21^ is of the form 

fe{p,q)=(^eip)Peiq), (27) 



14 



We remark that Euler's number e in (I25p is inherited from the iV! term in (|20p . 
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where Pe{q) solves the following fixed point equation on q space, 

^ ^ exp {-^,{p,)-^ Ujyjq, q)ps{q)d^q) ^^g) 
/a^^P {-^eipey^ J^UA{q,q)peiq)d^q) dq 

with "dsip) given by 

IMp) =e-jj lU^iq, q)p{q)p{q)d'qd% (29) 
and where cre(p) = a{p£){p), with a{p){p) defined whenever {}e{p) > 0, by 

a{p){p) = (27r^,(p))-iexp(-i|p|V^9.(p)). (30) 

Evidently, every minimizer of ^b(/) over 2te factors into a product of a 
Maxwellian on p space and a purely space-dependent "self-consistent Boltz- 
mann /actor. 'o However, the Maxwellian in ( 127|) is not autonomous from the 
Boltzmann factor in (127|) . as is manifest by the functional dependence of the 
(rescaled) temperature "i? = i^eiPe) on p^, see (130!) . For a subset of e values the 
minimizer of ^b(/) over 21^ may not be unique, but all minimizers produce the 
same asymptotic formula ( 123|) . In such a case of non-uniqueness of minimizers, 
they always seem to constitute either a finite set (typically a first order phase 
transition) or a continuous group orbit of a compact group (e.g., when A is in- 
variant under 5*0(2) or 50(3) and a minimizer breaks that symmetry), to the 
best of our knowledge; this seems to cover all physically relevant possibihties. 

In addition to the minimizers of ^b(/) there may be non-minimizing critical 
points of ^b(/) satisfying f l7r|) -fl5Ul). but these are irrelevant for fl2^ . 

Our Theorem [31 proved in section 5.3 with input from section 5.2, charac- 
terizes the Vlasov limit — )■ oo of the marginal measures 

>f ^(d'^'^X) = /if ^(d'^'X X (R^xA)^-"), n = 1,2, ... {n fixed) (31) 

in terms of the /g. We note that the object of interest in (mathematical) 
physics is not ([2]) itself but only the collection of its first few marginal measures 
(PTj) . To state our theorem, we introduce *P'^((M^ x A)^), the permutation- 
symmetric probability measures on the set of infinite sequences in M'^ x A. A 
theorem of de Finetti |deF37] . Dynkin |Dyn53] , and Hewitt-Savage |HeSa55] 
(see also |E1185j . App.A.9.) states that ^*((M^xA)^) is uniquely presentable 
as an average of infinite product measures; i.e., for each p, e ^^((M^xA)^) 
there exists a unique probability measure z/((ir|/i) on *!p(M^ x A), such that 

>(d3>d3"g) = [ r®"(d2pid^gi • • • d^p„d^g„) u{dT\p) Vn E N, (32) 



^^The expression conventionally known as "Boltzmann factor" results when Wa = so 
that Ua{<1, q) = V\{q) + Va{ci), i.e. for the perfect gas acted on by an external potential Va. 
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where "ju is the n-th marginal measure of /i, and r®"(d^pid^gi ■ ■ ■ d^p„d^g„) = 
r(d^Pid^gi) (g) ■ ■ • ® r(d'^p„d^g„). Equation (132|) is also the extremal decompo- 
sition for the convex set ^'*((M^ x A)^), see |HeSa55] . 

Theorem 3. Under the same assumptions as in Theorem\^ consider ^ with 
Hamiltonian ( fTPj) as extended to a probability on (M'^xA)^. Then the sequence 

{/^Tvaei^GN tight, SO onc can extract a subsequence {lJ'^^2^^^}Nm such that 
lim >5fJ^lV^>d=^"?) = >e(d^>d3"9) e ^'((M^xA)") Vn G N. (33) 

r/ie decomposition measure v{dT\fiir) of each such limit point fie is supported by 
the subset ci/^(]R'^xA) which consists of the probability measures Ti;{d^pd^q) = 
f^[p,q)d^pd'^q which minimize the H functional ^b(/) over Sl^. 

5 Proofs 

We have stated our Theorems 1,2,3 entirely in terms of the familiar quanti- 
ties of kinetic theory. These are the one-body density function fs{p, q) which 
minimizes Boltzmann's if-function ^(/) under the familiar energy functional 
constraint £ (/) = e, and the system-specific Boltzmann entropy per particle 
Sj^{e) which is given as the negative of Boltzmann's if-function evaluated with 
/g. However, in this format our theorems give essentially symmetric weight to 
the p and q variables, which ignores the fact that the p-space integrations in- 
volved in (I3T]) and ( 120|) can be carried out explicitly in the same fashion as for 
the perfect gas. As a consequence the problem reduces to studying the large N 
asymptotics of the expressions which result from these p-space integrations!^ 
In fact, all the hard analytical work goes into controlling the q-space integra- 
tions. This is certainly the case as far as the entropy per particle goes, yet also 
each minimizer of ^b(/) over the set 21^ is uniquely determined by p^, which 
signals that all of our Theorems 1 to 3 will be essentially straightforward corol- 
laries of theorems about certain qr-space expressions. Those theorems take a 
less familiar form, presumably, which is why their statements have been rele- 
gated into this section where we prove Theorems 1 to 3. 

5.1 Proof of Theorem [1] and its two upgrades 

To prove Theorem [1] we first formulate and then prove an upgraded version 
(Theorem [TJ'"), whose proof also proves Theorem [TJ 

^^AU Boltzmann needed for this was that (1 + x/n)" x e^; cf. |Bol96j . part II, ch. 3. 
Of course, things are not quite as straightforward with an irreducible W\ ^ 0, or else 
Boltzmann would not have had to have W\ ^ excluded from his analysis. 
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5.1.1 Theorem [T}^ and its proof 

Carrying out tlie p integration£l| in il' (^v) (£) given by iD, with h[^^ given m 
(fT9|l . Boltzmann's ergodic ensemble entropy (!20|) becomes 

V-,(£) =ln((HZ^|§3A^-i|v^' (£)) (34) 
with [S'^^^^l the standard measure of the unit 3A^ — 1 sphere S^^"-*^, and with 

= w^. y(£-/r(gi,...,q;v)) ^ xs^ji^^^^,f\ (35) 

where we introduced the interaction Hamiltonian 

/r(qi,...,q^) = 5^5^f/A(q„q,). (36) 

l<i<j<N 

Implementing the Vlasov limit scaling, i.e. setting £ = N'^e with e > Eg > 
0, recalling that IS^^^-*^] = n^^^"^ /T{3N/2), and using Stirling's formula for 
Euler's F function, we obtain the following asymptotic expansion for (|34l) . 

S^(N){N^e) = -iVlniV + A^ln (^|A| (^e)^^^) +0(lnAr) 

+ In j{l - ^lf\q„...,q^))^ ' A(d^^g). (37) 

where (■■■)+ means the positive part of (■■■)) moreover, A(d3^g) is the A^-fold 
product of the normalized Lebesgue measure X{d'^q) = \A\~^d^q on A. For 
brevity we wrote |A| for the volume Vol (A) of A. 

When = in A^, then becomes the Hamiltonian of the perfect 
gas without external fieldsQ abbreviated as K^^^ (for kinetic Hamiltonian). 
In this case the second line in fl37|) vanishes, and fl37|) becomes the asymptotic 
expansion of the entropy of the spatially uniformly distributed perfect gas, viz. 

S^(N){N^e) = -iV In iV + iV In (^|A| (^e)^^^) +0(lnAr). (38) 

The coefficient of the 0{N) term in (138|) gives the system- specific Boltzmann 
entropy per particle of the spatially uniform perfect gas, which we denote by 

s,Ae) = ln{\A\{'-fef'). (39) 



^^It is understood that d^'^ X etc. now involves the p variables used in 

Led 

configurational integrations in (|37p are still over 



^^It is tacitly understood that the cutoff provided by remains effective, so that the 
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Whenever interactions ^ of the admitted type are present, Theorem [T] 
follows if we can show that the second line in (|3711 is 0{N) and so contributes 
additively to the system-specific Boltzmann entropy per particle, and provided 
it has the right monotonicity and regularity. This is expressed in 

Proposition 1. Under the assumptions stated in TheoremUl there holds 

Jim i In y(l - ^/f " \{^^\) = s,,{e). (40) 

The function e i— )■ s^^j{e) is continuous and increasing for e > Sg > 0. 

This concludes the pretext for our first upgrade of Theorem [H stated next. 
Theorem 1.+ TheoremUl holds, with 

where s^^^^{e) is given in ^3^, and s^^j{e) in (J^- 

Proof of Theorem : 

Clearly, Proposition [1] and formula ( |3711 imply Theorem [1] and the splitting 
of the system-specific Boltzmann entropy per particle {e) in fl22l) into a sum 
of a kinetic and an interaction component, fHT]) . Proposition[T]also adds a piece 
of information about Sj^^j{e) which does not just re-express what is stated in 
Theorem [H In fact, by the known strict increase of e i-^ Ine, the increase 
oi e ^ Sj^^j{e) implies the strict increase of e h> ^^(e), but the increase of 
e I— Sj^^j{e) does not follow from the properties of e ^ \ne and the strict 
increase of e i— )■ s^(e). So Theorem [1}^ holds and extends Theorem [H □ 

Proof of Proposition \^ 

By hypothesis (-^6), [/a is bounded continuous on Ax A, so we can write 

iV-2/f)(gi,...,g^) = ^if/A(q,g)Ag,,(g)A«,,(g)d3gd3g 

IW^{q,q)A%{q)d\ (42) 

and we may abbreviate the first term of r.h.s.( H2l) in bilinear form notation, 

II lUj,{q, q)A%iq)^%^m'qd^ ^ {A%, , A%,). (43) 

The above integrals extend over M'^, and we set I^\qi, ...,q^) = cx3 as well 
as (A^Ijv) 5 A^l^)) = oo if any ^ A. Also by {HQ), the term in the second 
line of r.h.s.(H2|) is 0{N~^) for all X^^^ G A^. Recalling our claim (which we 
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promised to prove) that the hmit — > oo for the ensemble does not change if 
the Hamiltonian is changed by an additive term of order 0{N~^) relative to 
the leading terms, we now introduce the configurational integral 



Tf)(5) = ln/(1 



A 



(1) 



A(d=^^g) 



(44) 



for all > Nu{e) (to be defined). Note that the integral (14^ is gener- 
ally not well-defined for all G N because (A^l^, ^A^Iat)) is bigger than 
^_2jW^g^^ ^j^g absolute value of the second line of r.h.s. p2|) . which 

reads precisely A^~^ X^^i ^U^^qk, Qk)- while this term = 0(1/A^), when 
is not large enough then it is possible that ( A^^j^) ? A^(jv) ) > ^ everywhere 

in A^, in which case the integral in ( l44l) vanishes, and its logarithm = —oo, 
then. Yet, when A^ > Nu{e) the integral (jHj) is well-defined, and we conclude 
that (modulo the proof of precisely the just re-uttered claim that 0(1/A^) con- 
tributions to the Hamiltonian drop out when A^ — )■ oo) our proposition [1] is 
proved if we can prove the following proposition. 



Proposition 2. Under the hypotheses on U\ in Thm. when Nu{e) then 
N',,(c) + o(N) . (45) 



The function e i— )■ 7a(£^) is continuous and increasing for e > Sg > 



Proof of Proposition \^ 

We will establish uniform bounds and super-additivity estimates. 

For < n < A^, we set X(^) = (X("), F^^"")), which also defines 
We note the convex linear decomposition 



f)ASiv-„)(9)- 



(46) 



Since ?7a > is the kernel of a bilinear form which is positive definite when 
restricted to the set of probability measures on A, Jensen's inequality gives us 



We of course also have 1 



N ' 



±±.y{]V-n) ) £±y(JV-n) / • V^' ) 



^ + (1 — ;^), and so we conclude that 



e\^XiN) ^ i^X(N) / ) - 



N 



l_i/A« A^') \ 



1-KA?L^),aS--")> ) • (48) 



Next we recall that, if is some function on a domain and if S((y9+) 
denotes the support of its positive part, and is the characteristic function 
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of S ((/)-(.), then the inclusion S(<y9+) fl C T,{{ip + for any two such 

functions (f and ^ yields the estimate 



E{v+)AE{i)^) 



(</'++^+)x.(.+)X,(,^)- (49) 



Qpt w - 1 - i/A^^^ A^'^ \ 
aei - ^ i 



.(1) 



and 



e\— y(iv-") '— y(iv-")/ • 
ying next the classical in- 



Then inequahty ( 1491) applies to r.h.s. ( H8l) . App 
equality between the arithmetic and the geometric means of any two positive 
numbers A and B, viz. a A + (1 — a)B > A°'B^^~°'^ for any a G [0, 1], we get 



1 - i/A^'^ A^'^ 



(50) 



We now use (150]) to estimate r.h.s. (jS]). For this, let > A^c/(^) and let 
Nu{e) < n < N — Nu{e). Noting that the resulting integral over factors 
into two integrals, one over A" and another over A^~", and working out the 
powers, we find 



lny(l-l(AS.),AS.,>); A(d^-g)> 



(51) 



ln/(l-i(A^|„,,A^|„>)^ ^^A(d^) + 



3n _ n. 

2 AT 



.(1) 



3(N-n) ^ ^ n_ 



N 



A(d3(^-")g), 



where we also relabeled the integration variables under the second integral on 
r.h.s. (15T]) from yC^-") to X^^~"'\ Noting next that < < 1, we resort again 
to Jensen's inequality, this time w.r.t. the A measures in the two integrals on 
r.h.s. (IST] ). Also using ln(- ■■)"' = aln(- ■ ■ ), we arrive at 



ln/(l-l(AeLn,A«,,>) ' 'Xid'\)> 



3N 



■1 



1 + 



2~2n/N\ i„ If, 1/^(1) ^(1) 



(52) 



3n-2 J 
2n/N 



In 1-K4S„,4S„> A(d^) + 



' 3(Af-n)-2^ 

Formula (152|) writes shorter thusly, 

Tr(-) > (1+ ^) Tr(-) + (i + (53) 
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(e) is almost super- additive. 

To be able to create a properly super-additive function we establish upper 
and lower bounds of £ i— )■ T^l\e) which are linear in i, whenever i > Nu{e)] we 
will need those bounds with i G {n, N — n}, with i > 1. As a by-product, the 
upper bound with i = N will also guarantee convergence of the constructed 
super-additive function. 

The upper bound is trivial. Recall that by hypothesis ( A^^^, , A^|^j ) — 
for all i eN. So for i > Nu{e) and e > Eg > we find 

ln/(l-H4S,,AS.>)j"A(d^^g)<0. (54) 



3e-2 



As for the lower bound, we distinguish two cases, (a): (lAI^-*^, lAI^-*^) < e, 
and (b): (lAI^"*^, |A|^^) > e. In case (a) we apply Jensen's inequality w.r.t. A 
to the convex map x {1 — a;)^ (for 9 > 1), and also use ^ — 1 < ^, to get 



In 
In 



1-1(aS.,aS,>V A(d^) 



2 

3<-2 



> 



Ve J lUA{q,q)X{d'q) - jj i^A(g, g)A(d^g)A(d= 



(55) 



and r.h.s.(l55l)> —C > — oo when i > icriti^) (given Ua), with C > inde- 
pendent of i. Since the interaction entropy exists when i > Nu{e), clearly 
^crit > Nu{e), but after at most an adjustment of C, we can conclude that 
l.h.s.(l55])> -C > -oo when £ > Nu{e), with C > independent of i. In 
case (b), inequality ( j55l) is still true but now trivial, for r.h.s.(l55l)= — oo 
for all £ > 1, then. So instead we now proceed as follows. By hypoth- 
esis (-f^6), the bilinear form ( A^j^, , A^|^j ) takes its minimum £*g{i) > £g- 

Clearly, e*(£) = eg{i) + 0{i~^), where eg{i) := mmi-^li!;\q^, q^), and since 
^g{^) < £g (as proved in Appendix A), we have that e^ii) < + 0{i~^); of 
course, we also assume that i > Nu{e) so that £*g{(^) < £■ By permutation sym- 
metry there are many equivalent minimizers, but possibly also several distinct 
permutation group orbits of minimizers. We pick any particular minimizer Qg 
and let q^g \ G A denote the k-th coordinate vector in Qg^"*. By {H6) again, we 
can vary all the in the minimizing configuration a little bit, say, each 
in Bs{q^g\) fl A, where Bs{q) is a ball centered at q, with radius 6 > inde- 
pendent of k and i but chosen small enough (given e) so that ( A^j^, , A^j^, ) 
does not change by more than {e — eg + 0{i~^))/2. For brevity we write Bs[k] 
for Bsiql%); 

XBs[k] be the characteristic function of Bs[k]. We use that 
X{d^qk) = XBs[k]X{d^qk) + XBi[k]X{d^qk) where Bl[k] = A\Bs[k] is the com- 
plement in A of -B<5[fc], then use that both terms in this decomposition are 
non-negative so that we get an upper estimate by dropping the contribution 
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from XB^[k]^{d^<lk) for each k. After this step the restriction to the positive 
part of (! — -(., . )) is eventually tautological when i is sufficiently large so 
that the 0{i~^) term has gotten sufficiently small. We next apply Jensen's in- 
equality w.r.t. the probability measure Yli<k<eilBslk]nA K(^^Q))~^XBs[k]K(^^Qk) 
to the convex map x+ i-)- (for > 1), finally recall that < Eg < e, and get 



(1) \ 

■xw I 
(1) a(i) 



2 

3«-2 



> 



1 - i/A^'^ A^'^ \ 



A(d3^g) 

n XB,\k\K^qk) 



\<k<l 



2 

3f-2 



> 



iairf-/(i-KAS„.AS„>) n Eisf^w^d'* 



> 



i<k<e 



\Cs\^' (1 - 1(1 + ^) + o(r^)) > c > 

for i large enough; here 



mm 



s(<7')nA 



X{d^q) > 0. 



(56) 



(57) 



In summary, our list of inequalities (154|) . (!55|) and (!56|) . and the finiteness 
of the number of i until "£ is large enough," establishes that when £ > Nu{e), 
then for some ^-independent constant C^, > 0, 

-(f -l)a<Tf(5)<0; (58) 

incidentally, our ( l55l) and ( l56l) produce an upper estimate for Nu{e). 

Recall that in this proof we assume that ^ Nu{e), and that Nif{e) < 
n < N — Nu{e). With the help of f l58|) . for i e {n, N — n}, we conclude from 
fl53|) that there exists a C G M independent of n and such that 



T 



.)>Tf (5) 



(59) 



Adding that constant C to both sides of the inequality (l59ll shows that A^ i— )■ 
Tj^''(e) + C is a super-additive function for all e > > 0. And using ( l54l) 
with A^, we also see that Ar-i(Tf^ (e) + C) is bounded above, and so, by 

standard facts about super-additive functions, (T^^\e) +C) converges as 
A^ oo, 



lim N-^ 

AT-j-cxD 



and since A^ 



: supAr-ifTf)(£) + c), (60) 
0, we conclude that A^~^T^^(e) converges as well, i.e. 



{N) 



(61) 
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This proves (H5|l . 

To prove continuity of 7^(^)5 we establish upper and lower bounds on the 
derivative of the functions e 1— N~^T^^\e) which are uniform in > Nu{s). 
Differentiating the functions e 



-iTX'^(e) + (I - ^) Ine, we obtain 



(3 _ n i 



/(I 



i/a(i) a(i) \y 



'A(d3^g) 



_/(i-i(AS.,,A(;u»7-^(d3^) 



(62) 



To get a lower bound, we split off a factor (l — -^(A^^jv) ) A|^|jv)))+ in the 
integrand of the denominator of r.h.s.(l62l). and using that s > Eg > the 
positivity of the bilinear form now gives (l 

iTrV)>(f-^)^[i-i] 



0; 



^xiN) , Ay^)))^ < 1, and so 

(63) 



incidentally, this shows once again monotonicity \ oi e ^ 7a (^)- To get an 
A^-independent upper bound to ([62]), note that ^-2 = (^-l) (l- -^^) 
and that < (l — ^^^2) < 1 for > 1, then apply Jensen's inequality w.r.t. 
A to pull the power (l — out of the integral in the numerator, then note 

a cancellation versus the denominator. Since < (l — ^^_2 ) < 1 for > 1, 

(64) 



1 



e) < 



TV 



3iV 1 


2 -1 

3JV-2 


); "A(d3^) 


- 1 



whenever A^ > Njj{e) (so that the integral is non-zero). By the first inequality 
with i = N, the r.h.s.(1M|) is bounded above independently of A^. The 



m 



continuity of e H- 7^(5) follows. 

Proposition [2] is proved. □ 
To complete the proof of Proposition [1] we still need to show that the 
omission of /|f/A(q', g)A|^|jv) (q)d^g from P2|) was justified. This is now 
straightforward. By hypothesis {H6), U\{> 0) is a bounded continuous func- 
tion on A X A. So there exists an A^-independent constant -B > such that 



0< I 'Mq,q)A%,{q)d'q<B 



(65) 



as long as X^^^ E . Thus, and abbreviating the expression in the second 

(66) 



line on r.h.s. (1571) by S (N)(N'^e), we have the two-sided estimate 



But 



e+BN-^ 



Tf^'(0 



d<; < BC, 



(67) 
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the last inequality by and by the first inequality in (158|) . with i = N, and 
by e < <; < 25. So we conclude that for any i? > we have 



hm ±T'^^{e + BN-')=^,is). (6^ 
Hence, and by 



hm ±S^,.,{e)=^,{e), (69) 



and Proposition d] is proved, with s^^^{e) = 7a(£). □ 
This also completes the proof of Theorem [H □ 

5.1.2 Theorem [1}^"'" and its proof 

Ruelle's proof |Rue69] of the traditional thermodynamic limit for f l20|) per 
volumq^ proceeded along somewhat different lines, and when adapted to the 
Vlasov scaling it yields an interesting alternate proof of Theorem [T] which 
characterizes ^^(e) in terms of a variational principle (VP) involving s^^^{e) 
and yet another (auxiliary) "interaction entropy," which we denote by ^^^^-(e). 
For technical reasons we now need to assume that £g > (rather than Eg > 0). 
So, following Ruelle |Rue69j we introduce the configurational integral^ 



S/f)(£) = yX|,(iv,<,|A(d^^g). (70) 

Up to a purely numerical factor, (1701) is quasi the "3A^/2-th derivative" w.r.t. 
£ of \l',(]v)(£), the first derivative of which is given in (|35|) . For convenience 

A 

we rewrite (ITOl) . with £ = N'^e, as 

X n 

(71) 



Proposition 3. Assume the hypotheses of Theorem but now let Sg > 0. 

Then the following limit exists, 

hm ^lnH^,.)(iV2e)=^^,^(e), (72) 

N^oo -"A 

ands^^j{e) < is an increasing, right- continuous, function of e > Eg. 



Actually, Ruelle discussed the entropy of a regularized microcanonical ensemble measure 
|Rue69| . In |Kie09a| the author showed that a minor modification of Ruelle's approach 
establishes the thermodynamic limit for (|20p per volume without regularization. 

^°Instead of the normalized Lebesgue measure A(d^^q), Ruelle |Rue69| uses Nl^^d^^q 
which gives equivalent results in the thermodynamic limit; not so in the Vlasov limit. 
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Proof of Proposition 

Simplest things first, we note that S {jv)(£) < 1 (obviously), which proves 

-"a 

that In S (jv) (A^^e) < for all A^, and so s^,j(£) < whenever this limit 
exists. The proof that this limit exists and is a monotonically increasing right- 
continuous function oi e > Eg > consists of two main steps. 

First, as in our proof of ThmJH we temporarily replace N~'^I^\q^, q^) 

by ( A^Ijv) ; Ax(jv) ) ill dSl) ciiid study its logarithm. For this we need once again 
to assume that N S> Nu{e). Inspection of our proof of Proposition [2] reveals 
that we can recycle inequality fl50|) . take its vanishing power, integrate and 
take logarithms, and for Nu{e) < n < N — Nu{e), in place of flSTl) we now find 

1^ /(^ - i^xlm ,A_%))\{d'\) > (73) 
ln|(.-(A«,,A«,>)V"g) + 

which proves super- additivity of N^\nJ{e - (A^J^, , A^|^)))° A(d3^g) with- 
out further ado. Furthermore, since (■■■)+ is either 1 or 0, we conclude that 
In /(e — (A^Ijv) , A^|jv)))° '^(d^''^9) < 0. This upper bound and super- additivity 
now yield that the following limit exists, 

jim iln |(. - {A%, ,A%))\{d'^q) = S.Je); (74) 

moreover, s^,/!^) ^ is monotonic increasing, since l.h.s.f l74p is. 

Next we would like to prove continuity of s^^j{e) as function of e and then 
conclude the proof as at the end of the proof of Theorem [T|, but so far a proof 
of continuity of eluded us. Fortunately we can bypass this obstacle 

because s^^j{e) is a monotonic increasing function of e. We define 

Sa,7(^^) = inf Sa,/!^;^) (75) 

X>1 

and show that 

jim iln - ^/r(9i,...,9^))° A(d=^^g) = S,,(e+), (76) 

which proves Proposition [31 with 5^,^(5) = s^^j{e~^). 

To accomplish this, we recall P2|) and P3|) and rewrite f l7T|) as 

^.iMN'e) = J{s- {A% , A%) + j^{A%,))\{d'\), (77) 
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where we also introduced the abbreviation 



(Ag,,> = I \U^{q, q)A^X{q)d\ (78) 
Since now £g > 0, there exist constants B_,B satisfying 0<B_<B <oo so that 



B<{A%)<B. (79) 
But then, for all > Nu{e) big enough, we have 

^ In {N'e) > {e + N-'B) + o(l) > + o(l) (80) 

where o(l) — ?■ as — > oo. So 

hm inf ^ In S^,^) (N^e) > ~s,, {e+) . (81) 

On the other hand, for all > Nu{e) we also have that 

ilnS^(^,(Ar2£) < + AT-i^) +0(1), (82) 

and so 

lim sup ^ In S (^) (AT^e) < 5^^^ (5+) . (83) 

The estimates ( 18T|) and (183|) prove (1761) . 

So s^,j(£:) = s^^j{e~^). Of course, s^,j(£^) = 5a,/ (^) all £ which are points of 
continuity of Sj^^j{e), and the two functions share their points of discontinuity. 
At such points s^,j(£^) is right-continuous and may or may not agree with 

Proposition [3] is proved. □ 
We are now ready to state our second upgrade of our Theorem [H 

Theorem 1.^^ Under the hypotheses of Proposition\^ TheoremUl holds and 
the system- specific Boltzmann entropy per particle {e) given in 1^2^ satisfies 
the variational principle 

s^{e)= sup (s^,^(xe) + s^,,([l -x]£)). (84) 

0<a:<l ^ ^ 

Proof of Theorem : 

Integration by parts yields, for any £ > and e > Eg, 

1 rW^S^^37V^_ r [(. 1 







1 - ^ ) , A(d-g) = / ^^IT^ ] . A(d-g)dx^ (85) 
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where we suppressed the arguments {qi, q^^) from 1^ {Qi, Qn)- Setting 
£ = ^ - 1, recalhng ([71]), and using that A^"^ In - l) ^ 0, we find 

s^,(e) = hm ^In / EAN){N^[l-x]e) x'^-^dx. (86) 
Jo ^ 

Proposition [3] and Laplace's method (cf. sect. II. 7 in [E1185] ) now yield 

SA,ii^)= sup (|lna; + s^,,([l -x]^)); (87) 

0<x<l ^ ^ 

note that ( 187|) implies that e s^^^^e) is continuous even when s^,^(e) is not. 
Recalling next the definition (159]) of ^^(e) as well as ( 1^ of Theorem [TJ*', we 
see that Theorem [T]^"'" is proved. □ 
We end this subsection by pointing out that our method of proving The- 
orem [U^"'" not only avoids the regularization of Dirac's 6 measure, we also 
tackled the map £ t-)- S{S,) directly rather than its inverse S i— )■ -E(S) |Rue69] . 
The strategy to tackle § ^ E{S) is due to Griffiths [Gri65j . 



5.2 Proof of Theorem 



Since formula (1371) holds also under the assumptions {H1)-{H5) on the inter- 
actions, and since it is well-known that the system-specific Boltzmann entropy 
per particle of the perfect gas ( 139]) minimizes Boltzmann's H functional under 
the constraint of prescribing the value of the kinetic Hamiltonian, it suffices 
to study the interaction entropy of Boltzmann's ergodic ensemble, 

S^M^) = \nj {1 - llf\q,,...,qj,))f-'X{d'\). (88) 

Note that ( 188]) is non-positive, and under hypotheses {H1)-{H5) we also have 



1- ^1^/^(91, ••■,9^))--^A(d^^g) 



> 



l<k<N 



> 



ICJs^w / 1 



(1-^1 + ?))>» 



l5 •••5 



n 

l<fc<Af 



/si-[fc]nA^(d^9) 



^9) 



where again Cs is given in (|57]) . but now with 6{6) independent of k and 
chosen so that N~'^I^^\q-^^, q,^) < eg{N) + {e — eg)/2 when the q^. 
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vary in Bs^qg jS) fl A, where (q^^i, ...,9^ ^v) is a ground state configuration for 
I^\q-^, q^) witli a fat neigfiborhood, whicfi exists by (iJ2)&(if3). We also 
used tliat £g{N) = min N~'^I^j^\qi, ...,qjv) — (see Appendix A). So 



2 



3Af-2 



^,(iv,(£) > In (iC^li (1 - i (1 + ^))) > -oo (90) 



for all > 1. The estimate fl90|) guarantees the existence of limit points of 
the (negative) interaction entropy per particle as A^ — )■ cxd. We want to show 
that the interaction entropy per particle actually has a limit and characterize 
the limit by the variational principle stated in Theorem O 

We begin by characterizing ( IHHi) by its own maximum entropy principle. 
We introduce the quasi-interaction energy of g^^^ E ^*(A^), defined by 



(91) 



whenever supp f?*^^^ C supp [e — N ^/i^'*)^; else we set Q,^^^^^ (f?*-^-*) = —oo. 
The entropy of g^^^ relative to gip^ E ^''(A^) is defined as usua|^ by 

^^^^ i^n^iV) = - /in (^^) ^^(^Hd^^) (92) 



if is absolutely continuous w.r.t. the a-priori measure gLp \ and provided 
the integral in fl92|) exists. In all other cases, ^L^^-* (f?^^''|f?ap^'') = — oo. Finally, 
we define what we call the interaction entropy of g^^^ by 

^ ^(N) + (^(A.)) . (93) 

We are now ready to state our variational principle. 

Proposition 4. For e > Eg > 0, the interaction entropy functional (E 
achieves its supremum. The maximizer is the unique probability measure 



(i-^/r(g„...,g,))¥-d3-g 



Q^^lid'^'q) = -^^ E (r ni:°°)(A^); (94) 



thus 



Moreover, 

s'j/X&l) = SMN'e). (96) 



'^Our physicists' sign convention of relative entropy is opposite to the probabiHsts' one. 
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Proof of Proposition 

Under our hypotheses on I^'' the measure g^^2 is absolutely continuous 
w.r.t. A and bounded whenever e > Sg, so the standard convexity argu- 
ment due to Boltzmann |Bol96j . cf. |Rue69l IE1185] . applies and shows that 
SjJ^{g^^^) — SjJ^{g^^2^) < 0, with equality holding if and only if g^'^^ = g^^2^. 
Identity (l96l) is verified by explicit calculation. □ 

Since ultimately we are interested in the limit — t- oo of our finite-A^ 
results, we recall the formalism of probabilities on infinite sequences A^, as 
encountered already in section 3 for (A x M^)^. Thus, by ^*(A^) we de- 
note the permutation-symmetric probability measures on the set of infinite 
exchangeable sequences in A. Let {"gjneN denote the sequence of marginals of 
any g € <P"(A^). The de Finetti |deF37j - Dynkin |Dyn53| - Hewitt-Savage 
|HeSa55j decomposition theorem for *P'^(A^) states that every g G ^'^(A^) 
is uniquely presentable as a linear convex superposition of infinite product 
measures, i.e., for each g G ^''(A^) there exists a unique probability measure 
<;{dp\g) on *P(A), such that for each n eN, 

>(d'"g) = / P'^^ld'gi ■ ■ ■ d'qn) ^{dp\g), (97) 

where '^g is the n-th marginal measure of g, and where p®"'(d^gi ■ ■ ■ d^g„) = 
p(d^gi) X ■ ■ ■ X p(d^g„). Also, fl^ expresses the extreme point decomposition 
of the convex set ^'*(A^), see |HeSa55j . 

Next we would like to formulate the N = oo analogue of (p3l) . but the 
naive manipulation of the formulas is not recommended. The functional 0,^^^^ 
is well-defined by (!9T]) and its accompanying text for all G N; however, since 
our conditions on I^\q^, ...,q]^) allow it to be unbounded above when two 
positions q^. and qi approach each other (for example: Coulomb interactions), 
we find that Q,//^\p®") = — oo for all product measures p®"", but these are ex- 
actly the A^-point marginals of the extreme points of our set of exchangeable 
measures on the infinite Cartesian product A^. This obstacle can be circum- 
vented by noting that the finite- A^ quasi-interaction energy defined in fl9T|) and 
the line ensuing fl9Tl) is the monotone limit of a family of concave functionals 
in which the integrand function ln(l — x)^ (with InO = — oo understood) is 
replaced by ln(l — x)x{x<i-a} + [In a + (1 — a — x) / a]x{x>i-a}] thus 



"a^.x^ in = ^ /(in(i - i^/r)^{^r--(i-^)} 



where we omitted the argument {qi, q^) from for brevity, and 

aiX^ =limV,XH^^^"^)- (99) 
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We also define '^SjJ^{g^^'^) precisely like SjJ^{g^^'^) except that (l[^l{g^^'') 
is replaced by Q.^^/^ {g'-^^) . We have "^{^^^(p®") > -oo for all p e n 



£Mn£i)(A), and lim„;o "^J^ (p®") = whenever I/e ^ 1. By "^^^^ 



we 



denote the unique maximizer of °'SjJJ[g''^'^), easily proven to exist as done for 

Equally easily we find lim.,o "^j;,) ("^^.^J = S^J^^)- 
We are now ready to formulate the N = oo analogue of ( p3|) . 
To define the mean quasi-interaction energy of g E ^*(A^), we introduce 
the subset *^3^2(A^) C *P''(A^) for which the expected value of f/| is finite; 

i.e. J Ul{q, q')'^g{d^qd^q') < oo, where ^^(d^gd^g') is the second marginal 
measure of ^ G ^3^2 (A^). Also, by ^^72 (A) we denote the subset of q3(A) which 
consists of Lebesgue-absolutely continuous probability measures p for which 
/ U^{q, q')p^'^{d^qd^q') < 00, which implies {p , p) < 00; here we recycled the 
bilinear form notation ( H3|l for lower semi-continuous (rather than continuous) 
U\. If g E ^^2(A^), then the decomposition measure <^(dp|^) is concentrated 
on (A); this can be shown by adapting arguments from |HeSa55] : cf. also 
|MeSp82| . The mean quasi-interaction energy of g E *P^2(A^) is defined as 

Q, (f.) = hm hm '-"d^Pli-g). (100) 
We show that Q, ^^^(f?) is well-defined. By the linearity of "^g 1— )■ '^(iffj^i^g), 



the presentation (1971) yields 

= /""^Veipn^MQ), (101) 

and on (A^) the conventional law of large numbers for U statistics applies 
(see |Hoe48j ) and yields 

+ (lna + ^[l-a-7(P' P>])X{(,,,)>,(i_„)j 

Clearly, when a J, in f ll02p the "value" —00 is assigned to all p for which 
{p , p) > the a 10 limit is finite when {p , p) < e. We conclude with: 

Lemma 1. The mean quasi-interaction energy U00\) is well-defined and affine 
linear. For g G *P^,2(A^) having decomposition measure <;{dp\g) supported 
entirely by p for which {p , p) < e, we have UUU\) given by 



a,/,(f?)= / cij/M'^iMQ), (103) 
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where 

Q,,/,(p)^fln[l-i(p,p>]; (104) 

otherwise, Q, , (g) = — oo. 

The N = oo analogue of ( p2|) is the well-known mean (relative) entropy of 
g G *P*(A^), which is well-defined as limit 

^{g) = lim i^(")(>|A). (105) 

Here, 3?,^"^ ("f)|A), n G {0,1,...}, is the relative entropy of ""g, as defined in 
([92D; we also set ni^-^^ g\X) = for all k eN. The limit (fT05|) exists or is 
— oo. This is a consequence of the next lemma, which holds for g G *P'^(A^) 
or ^ G qp^(A^). If ^ = it is understood that A; < in gf\ 

Lemma 2. Relative entropy n i— 3?,*-"'' ("^|A) has the following properties: 

(A) Non-positivity.' For all n, 

r^(")("f?|A) <0; (106) 

(B) Monotonic decrease.' If n > m then 

(C) Strong sub-additivity.' For m, n < i, and k = i — m — n, 

^^'H'f?!^) <^^™n>l^)+^^"H>l^)+^^'H'f?l^)-^^"'H"'^l^)- (108) 

The proof of Lemma |5] is a straightforward adaptation from a proof by 
Robinson and Ruelle |RoRu67] (section 2, proof of proposition 1) for the 
standard-thermodynamic-limit problem to the Vlasov limit, studied here, cf. 
|Kie93j . 

The next lemma also has an elementary proof which likewise is an adaption 
from [RoRu67j . proof of their proposition 3, cf. |Kie93j . 

Lemma 3. The mean entropy functional UU5\) is affine linear. 

Lemma [3] in conjunction with the de Finetti |deF37] - Dynkin |Dyn53| 
- Hewitt-Savage |HeSa55j decomposition theorem for *P'*(A^) yields a key 
formula for the mean entropy which does not hold for the finite- A^ entropy. 
Namely, as a consequence of Lemma [31 the extremal decomposition of g yields 

^{g) = J ^^(p|A)^(dp|f?), (109) 

where we also set 3^(p|A) = 3^^^^(p|A). 

Lemma HJ also proved by adaption of a corresponding proof in jRoRu67j , 
proposition 4, ends the listing of properties of mean relative entropy (llOSp . 
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Lemma 4. The mean entropy functional is weakly upper semi- continuous. 



Finally we define the mean interaction entropy of g ^ *P'^(A^), 
By ffT09|) and ffT03D we have 

ll/ei0)= [ 5,,.(p)^(dp|^), (111) 

where we introduced the functional 

S^^^ip)^^ip\\) + (ll/e{p), (112) 

which is well-defined and finite whenever p G (^n£^ In 2,^)(A) and (^p , p"^ < e; 
else we have Sj^^{p) = — oo. Note that Sj^^{p) < 0, for si(^g\X) < and 
(li/e{p) < 0, the latter because J/a > by hypothesis. 

Because of (11 11 1) the problem of maximizing S_j^^{g) reduces to maximizing 
Sj^^ip) given in (11121) . 

Proposition 5. Sj^^{p) is weakly upper semi- continuous for e > 6g > and 

takes its finite non-positive maximum at a solution of the fixed point equation 

^ exp(-^-Hp)A^A(g,g)p(g)d^g) ^^^3^ 
J^exp {-^-^{p) J^UA{q,q)p{q)d^) dq 

where 

^9,(p) = |(l-l(p,p»5>0. (114) 



Proof of Proposition 

Since relative entropy 3?,(p|A) is weakly upper semi-continuous ( |ReSi80] . 
Suppl. to IV.5; |E1185] . chpt.VIII), and since the functional (li/e{p) is weakly 
upper semi-continuous as a consequence of hypothesis {H2) and the positivity 
of Ua, so is Sj^^{p). Since A is compact, Sj^^{p) now takes its maximum, which 
is non-positive because Sj^^^p) < 0, and finite (i.e. > — oo) because of the 
following. Let k pj-^^ in (*p fl (t^){A) be a minimizing sequence for {|p, p). 
Since e > Eg > 0, hj {H3) there is a K such that Eg < {p{k) , P(fc)) < ^ for all 
k>K. Then maxp5,/^(p) > 5,/^(P(i^)) = 3l(p(x)|A) + (li/e{p(K)) > -oo. 

Let q I— !■ Pe{q) denote any maximizer for Sj^^{p). Suppose (p^ , pe) > e. 
Then Cii/ei^Pe) = — oo, and because 3t(pe|A) < then also Sj^^{pe) = — oo. 
Therefore i^p^ , Pe) < e strictly, and since £ > 0, this proves (I114p . 

The standard variational argument now shows that the maximizer satisfies 
the Euler-Lagrange equation for Sj,{p), which is (11131) . □ 
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Corollary 1. The functional S_j^^{g) given in U10\) achieves its supremum. 
If Qe is a maximizer of s_j,^{g), then the support of its decomposit 



ion measure 

(,{dp\Qe) is the set of maximizers {pe} of the functional s ^^^{p) given in M12^) . 

Proof of Corollary \^ 

Abstractly, by Lemma H] and the linearity of the mean quasi-interaction 
energy functional, the mean interaction entropy functional S_j^^{q) given in 
f lllOP is weakly upper semi-continuous, and so achieves its supremum over the 
compact set of permutation symmetric probabilities ^^2(A^). 

Alternatively, by (11 111) and two obvious estimates, we have right away that 
■Sj/M = ^i/siP^) < sup5^/^(^) < max5,/^(p) = S^^^{pe), (115) 

so sViY>gS_ji^{Q) = max gS_j^^{g) = S_j^^{p^). Now let maximize S_j^^{g) and 
suppose that suppg(dp|p£) is not a subset of the maximizers {p^} of Sj^^{p). 
Then 

^i/eiSe) = / S {p) ^idp\ge) < max 5 (p) = 5^/,(p^), (116) 

SO g^ is not a maximizer — a contradiction to the supposition. □ 
We now relate the sequence of maximizers {^'j^^jArgN of {Sj^^} to the set of 

maximizers {p^} of s^^^. We begin with the maxima of SjJJ[g^^'^) and Sj^^{p). 
Proposition 6. We have 

^^j4Te{&)=^.jM- (117) 

Proof of Proposition 0.- 

For all a G (0, 1), we have 

"5;7;(«,Ki)>"5j;,Hpr)- (118) 

We compute 

-s\^^ApT) = N^^'\p.\X)+-ci\';i{pr). (119) 

Since (p^ , p^) < e, when a G (0, 1) is sufficiently small we have by (-^4) and 
Proposition [5] that 

J^Lr<^(/^n=|ln[l-KPe,Pe>]. (120) 



Hence, for all sufficiently small a G (0, 1 



30 



Thus 



iV^OO ' 



(122) 



for all sufficiently small a E (0, 1), and this yields the ffist desired estimate 

limjnfi5i;j(4^.i)>5,,,(p.). (123) 



Now consider ( jMll as extended to a probability on A . Since A is bounded, 
A is compact, and then the sequence {^'^l.jArgN is weakly compact, so 



lim '^gf}p = %e^^'{T) \fneN, 



(124) 



after extraction of a subsequence {^?*-^'^^''}AreN; note that the {"io^lneN form a 
compatible sequence of marginals. Furthermore, we have /^^^^^(d^g) = 0, or 
else Hi^Qe) = —oo, a contradiction; so "'g^ e ^^(A"). 

Following |MeSp82[ IKie93] we now use sub-additivity of relative entropy 
(property (C) in Lemma [2]) and then negativity of relative entropy (property 

(A) in Lemma [2]) (valid also with in place of g), and obtain 



< 



(125) 



where \_a/b\ is the integer part of a/b, and where m < n. Upper semi- 
continuity for the relative entropy gives 



limsup^(")("£[f])|A) < ^(")("^,|A), 



N^oo 



(126) 



while i 



— )■ ^. Hence, dividing f ll25p by N[N] and letting N ^ oo gives 

(127) 



lim sup 1^) < i^(")("^,|A) Wn e N, 

N-^oo 

and now taking the supremum over n (equivalently: the limit n — )■ oo) we get 
limsupi^(^)(4^,)|A) <^(^,). (128) 



Lastly, using ffT09D in ffT28D yields 

limsupi^(^)(4^,)|A) < / ^(p|A)^(dp|^,) 



N-^oo 



(129) 
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where <^(dp|^e) be the Hewitt-Savage decomposition measure for g^. For each 
p G supp^(dp|^£) we can choose a family of Q^^^p] E ^^{A^) satisfying 

hm "^(^Hpl =p®" (130) 
for each n G N, such that for each A^[A^], with G N, we have 

= / o'-''\p]<^{dp\gs). (131) 

In contrast to the de Finetti-Dynkin-Hewitt-Savage decomposition, this fi- 
nite decomposition is not unique, but this is immaterial. We remark 
that in the physically (presumably) most important situations, namely when 
supp (^(dpl^g) is either a finite set or a continuous group orbit of a compact 
group, then a decomposition fll3ip satisfying fll30p can easily be constructed 
explicitly, as shown in Appendix B. 

By (11311) . the linearity of the map g'^^'^ i— °Q, (f?*-^^) gives 

''^'f/Mi) = J''^TM%])^iMSe), (132) 
and by the concavity of the map / i— )■ 0.^1^/1 {g^^^) , Jensen's inequality gives 



lnri-ill(^)(f?(^))lx, , (133) 

+ flna + ih -«-il[(^)(f?(^))l)xr J, 

where 

U^'^Kg^'^^) = [l - N'') I |f/A(g,q)V^Hp](d^gd^g). (134) 
The weak lower semi-continuity of U/^ now gives 

hminf / |[/a V'^Wq > {P, P), (135) 
and since A^^^ — t- 0, we find for each convergent subsequence of measures that 

limsupi-diXH^^^^W) < I h[l - P>]x,,.,.K.a-.)> (1^6) 



N^oo 

+ (ln« + ^[l-«-KP'P>])x, 



(p ,p)>e(l-a)} 



for each a G (0,1). Now suppose that {p , p) > then r.h.s.( fT36l) .|, — oo as 
a J, 0, in which case by f ll36p and f ll32p also "Q-t/^ (f?7v2e) i as a | 0, 
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and by (!99|) and (!93|) and property (A) in Lemma [21 and then Proposition 
m this contradicts the lower bound ( pOl) . Therefore, {p , p) < e for every 
p G supp ^(dpl^e), and so, and recalhng (11040 . we conclude that 

limsup^cif/^{gfl)<l /ln[l-i(p,p>],(dp|&) = /(l,/,(p)^(dp|&). (137) 

The estimates (11371) and (I129p and two obvious estimates now give 

< limsupi^(-)(,!-)jA)+limsupiaiXH^!fi) 

N^oo N^oo N—^oD 

< / ^(p|A)?(dp|£-,) + / (li/e{p)<;{dp\Qe) 



^//e(p)^(dp|&) 

< max5//e(p), (138) 



and since this holds for each limit point we can drop the dot to get 



]imsnpj^sll\g^^l)<S,,M- (139) 

N-^oo 



By (11231) and (11391) . Proposition [6] is proved. □ 
By Propositions HI |5l and |6l the interaction entropy per particle for Boltz- 
mann's ergodic ensemble converges as follows, 

s,,{e) ^ hm j^S^^M^iN'e) = s.^Spe). (140) 

and s^,j{e) is characterized by its own variational principle expressed in Propo- 
sition [51 Moreover, by formula ([37]) . which holds under assumptions (-f^l)- 
{H5), the expansion f[23|) now follows, with 

5a(^) = 5a,a-(^) + 5a,.(^). (141) 

where s^,j^(£^) is given in f[39|) . and s^^j{e) in (I140p and Proposition [5] By 
Proposition [5l any maximizer p^ of Sj^^ (p) satisfies ( l28l) and ( l29i) . 

Lastly, one readily verifies that the just described s^{e) equals the negative 
minimum of Boltzmann's H functional over the set of trial densities 21^ = {/ G 
(*P n ii^ n £^ ln£^)(M^ X A) : e(/) = e}. This is done by explicitly carrying 
out the standard variational argument for H{f), taking the constraints into 
account with the help of Lagrange multipliers which are then eliminated with 
the help of the very functionals of pe displayed in Theorem [21 

This completes the proof of Theorem [21 □ 
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5.3 Proof of Theorem 3 



We begin with the observation that (!94l) is clearly the A^-th configurational 
marginal measure of ([2]), i.e. ( IMl) is (E]) with the Hamiltonian given by ( fT9l) . 
integrated over all the p variables in (|T9l) . Put differently, is the joint 
A^-point distribution on configuration space of an A^-body system with 
Hamiltonian f|T9|) chosen w.r.t. the a-priori measure ([2]) on (M^xA)^. Hence, 
our proof of Theorem [2] also proves the following weaker version of Theorem [31 

Theorem 3.^ Under the same assumptions as in Theorem 0, consider ^ 
for the Hamiltonian / fTPj) as extended to a probability on (M^xA)^. Then the 
sequence {Q^^2^}Nm of its configuration space marginals, obtained by integrat- 
ing over all the p variables in / fTPj) and given in (94\ ), is weakly compact in 

^*(A^), so one can extract a subsequence {f?^i^'^^'*}ArgN such that 



The decomposition measure <j(dp|^e) of each such limit point Qe is supported on 
the subset o/*P(A) which consists of the probability measures Pi;{q)d^q which 
maximize the functional Sj^^ (p) . 

Since each limit point of (p4l) is a convex linear superposition of infi- 
nite product measures on A^ consisting of "Boltzmann factors" Pe{Q) on A, 
satisfying (125]) with fl29p and maximizing the interaction entropy functional 
Sj^Xp)i since each such Boltzmann factor is associated with a unique 
"Maxwellian" (Je{p) on through fl30|) . each such Boltzmann factor thereby 
defines a unique Maxwell-Boltzmann distribution <ye{p)Pe{Q) on M'^x A given 
by the product of this Boltzmann factor with its associated Maxwellian. So 
the very decomposition measure ^(dpl^^) of each limit point Qs on A^ allows 
us to define a unique probability measure ps on (M^xA)^, viz. 



and this measure <^(dp|^e) on *P(A) can be mapped into a unique measure 
i/(dr|/ie) on ^(R^xA) which is concentrated on those r G *P(R^xA) which 
are of the form r(d^pd^g) = a{p){p)p{q)d^pd^q, with cr(p) given by fl5U]) and 
p satisfying (!28|) with (!29|) and maximizing Sj,{p), thus 




(142) 



in the sense that 




(143) 
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The corresponding infinite product measures ni<A;<oo '''{^^Pkd^lk) on (M^xA)^ 
are extreme points of ^''((M^xA)^), and so fll45p is tlie extremal representation 
of fis- So, having Theorem El" and its consequence ( 1145^ . all we need to do to 
finish the proof of Theorem [3] is to show that each such defined /ig is indeed a 
limit point of ([2]) under the stated hypotheses. 

To see this, we use that by Theorem [31" we already know that f ll43p holds, 

and we know the support of q{dp\g£). Writing ^g^^}^^^ explicitly gives 



= '""^ ^ . ' d'\ (146) 

where the integral in the numerator runs over the variables Qn+i to Qj^. For 
any and 1 < n < A^ we now write 

jr(g,,...,gJ = /^)(q,,...,qJ+/r'^Hq,,...,qJ+/f-")(g„^ (147) 

which defines /^"'^■'(q^, g^). Henceforth we omit the arguments from the 
Is to keep the formulas within sight; by fll47p the superscripts convey which 
variables are used. With the help of f ll47p we rewrite the integrands thusly@ 

/ An) , An\N) \ . . . 

= 1 ^ ^ ^ . Ul-4^^f'"^) -(148) 



el' 



Now (1 - (eN^y^I^^^)^ vanishes in an open neighborhood of configurations 
{q^, g„)oo for which I^^^ = oo, so that I^^ < oo on the support of fll48p . 
And for any configuration {q^, q^) for which < oo, we have N~'^I^^ — > 
as A^[A^] — )■ oo. Moreover, by our Theorem |3l" and its explication fll43p . 

we have that ^1 — -^^I^^^^^ ^-^d^'-^'^-'g, interpreted as a measure on the 

convex set of probability measures ^(A) with support in the set of empirical 
one-point "densities" with N — n atomso converges (up to normalization) to 



^^We are using that {fg)+ = f+g+ + f-g- for two arbitrary functions / and g, and that 
in our case / cannot be strictly negative if g is, giving (fg)+ = f+g+] in our case, / and g 
are the respective expressions between the two pairs of big parentheses at r.h.s. (ll48p . 

^•^If Ua is bounded continuous on A^, then we already know that we can rewrite Ij^^ as 
a sum of a bilinear and a linear form on *P(A) evaluated at a normalized empirical one- 
point "density" with N atoms; see (|42|) . If U\ is only lower semi-continuous this particular 
identification ceases to make sense, but happily we can always interpret as a linear 
form on the convex set of probability measures *P(A^) (cf. (|133p and (|135p ). evaluated at a 
normalized empirical two-point "density" with N atoms (O , and we note that any empirical 
two-point "density" with N atoms ([7|) is uniquely determined by its associated empirical 
one-point "density" with N atoms (jH]). 
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<^(dp|^£), and for any p in the support of <j(dp|^e) we have that N ^/|^ — )• 

(p,p) while iV-i/i"l^) ^ Ei<.<„JA^A(9fc,9)p(g)d^? when iV[iV] ^ oo. So 
for any p in the support of the decomposition measure <^(dp|^e) we have 

3JV 

/ j{n) j{n\N) / ] r \ 



l<fc<n 



with |'i?e(p) = £ — (p, p)- After this preparation, we now exphcitly compute 
the marginal "p!^'' and find 



(dv^) = ' ^ , — ^ — d^>d^"g, 



(150) 

where p„) = A^ X]i<A:<n ilPfcP- Using f ll47p we factor the inte- 

grands as in (I148p . though now we get 

' - 1 - tit.. 1 - i/f (151) 



and by following essentially verbatim the arguments which lead from 0148P to 
(11491) . we now find that for any p G supp ^(dp|^e), 

+ + ^ ^ -p. / ||p.P + A^A(q.,g)p(q )d^g 

(152) 

Our Theorem [3] is proved. □ 

6 Spin-ofFs of our results 

In this section we list a number of corollaries of our results. 

6.1 A weak law of large numbers / ergodic theorem 

Whenever ^b(/) has a unique minimizer over 2te, then necessarily all limit 



points in ( l33|) coincide, i.e. any fi^ = Pe- By the weak compactness of ^''(A ) 
(in product topology) we then in fact do have weak convergence, 

lim >5^2|^(d3>d3"g) = Xld^^pd^'"?) e ^'{{W^ x A)") Vn G N . (153) 
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Since in this case the decomposition measure u^drlfi^) is a singleton, the hmit 
/i£ = {"ju^jngN is of the form 

X(d3>d3"g)= H Mp,,q,)d'pA (154) 

l<A;<n 

with feip^q) = o'e{p)p£{q) as defined in Theorem [2l As discussed in |Spo91| , 
the factorization property (11541) is equivalent to a weak law of large numbers 
— or to an ergodic theorem, depending on ones point of view. Since the single 
particle momentum P and position Q of an individual A^-body system picked 
from Boltzmann's Ergode ([2]), with Hamiltonian (JT9l) . are random variables, 
any bounded continuous single-particle test function 6' on M^xA defines a new 
random variable © = 6'(P,Q), and so does its sample mean over a single 
A^-body system, 

1 ^ 

(®>^ = ivS^^^^'^^-^- ^^^^^ 

Theorem [3] in the special case (I154p implies that, for all such 6*, 

hm (0)^ = I e{p, q)ap, q)d'pd'q , (156) 
in probability. The generalization to n-body test functions holds as well. 



6.2 The Vlasov limit for other thermodynamic potentials 

A second corollary, or actually a whole family of corollaries, is the existence 
of the Vlasov limit for the thermodynamic potentials of the canonical and 
grandcanonical ensembles under the same hypotheses. We only discuss the 
Vlasov limit for the thermodynamic potential of the canonical ensemble. 

Thus, taking the Laplace transform of ([3]), i.e. multiplying by e~^^ and 
integrating over £, yields what is known as the canonical partition function, 

Z^,.m = ^ y"exp (-/3/7f )(XW)) d^^X (157) 

The Hamiltonian Hf\x^^^) is given in ( TT9l) . Clearly ( I157p factors as follows, 
Z^(^)(/3) = Zkw{I3)Z^,n){I3) (158) 

where 

Z^(.)(/3) = y"exp (-/3/f )(qi,...,q^))A(d=^^g) (159) 
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is the canonical configurational integral, with A(d^g) = |A| ^d^g the normahzed 
Lebesgue measure introduced in section [5l and A(d^^g) its A^-fold product, and 

ZKi^m = 4f /exp -.Pn)) d'^ (160) 

is the canonical partition function of a spatially uniform perfect gas in A, a 
Gaussian on the Cartesian product of the p spaces, which evaluates to 



(/?) = ^(2.tf)="'^ (161) 



here, we introduced N't} = f3^^, with ^ independent of A^, not to be confused 
with -i^e which is a functional of p. Since (3~^ receives the meaning of a tem- 
perature of a heat bath (up to the absorbed factor /cb), it needs to grow oc A^ 
to compensate for the growth of the system's energy £ oc A^^. Taking the 
logarithm of f ll57p gives what we call the canonical thermodynamic potential 
(canonical T-potential, for short )0 $ (jv)(/3). Using f ll58p and fll6ip as well 

A 

a.s (3 = ^ yields the asymptotic expansion 

$^(iv)(j^) = -N\nN + N\n {e\A\i27i^)^/^) + 0(\nN) 

+ lnZ,,.)(^). (162) 

Again, the A^ In A^ term is due to Gibbs' A^! and purely combinatorial in origin. 
In the absence of interactions (save the confinement to A) (11621) reduces to 

*xw(iv7^j) = -N\nN + N\n{e\A\{2n^f/^)+0i\nN), (163) 

the asymptotic expansion of the canonical T-potential of the spatially uniform 
perfect gas. The coefficient of the 0{N) term in fll63p is the system-specific 
Helmholtz T-potential per particle of the uniform perfect gas in A, denoted by 

0„,(^)=ln(e|A|(27r^)3/2). (164) 

The system-specific interaction Helmholtz T-potential per particle is defined 
by 

0,„(^)= hm ^lnZ^(.)(^). (165) 

The limit (11650 exists for Hamiltonians satisfying {H1)-{H5), as follows by 
corollary from Theorem [21 if {H2) is replaced by bounded continuity of the 
interaction, as explained earlier, then we can also infer the existence of the 



^^Multiplying the canonical T-potential by the temperature of the heat bath yields the 
negative of what is usually called the canonical free energy, which in the thermodynamic 
limit yields the Helmholtz free energy of the physical systems. 
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limit f ll65p from our Theorem [H The argument is quite standard, cf. |Rue69j . 
Namely, note that 

i In = F 1^ / e-^'+^('Me (166) 

where S{8,) is shorthand for Sj^{n){8.). Setting £ = N'^e and P = ^ and ex- 
panding S'(£) using (123|) (if [/a is bounded continuous on we can alternately 
use (1221) ), we find 



ilnZ^(.v)(^) = -lniV + ln 



^N{-d-^e+Sf,{e))+o{N)^^ 



^+0{^^). (167) 



Clearly, \\g\\^ — )■ \\g\\ao asA^— )■ oo, and so the following asymptotic expansion 
for the canonical T-potential results, 

'^<)(]^) = -N\^N + Nct>M + o{N) (168) 

with 

0j^9) = sup(-^-i£ + sj£)). (169) 

e>eg 

By (1I62D and ([l68]), we also have N-^ \nZ^i,N){^) 0a,/ (^)> with 

0a,.W=0a(^)-0a,kW- (170) 

This concludes our demonstration that the Vlasov limit for the system-specific 
Helmholtz T-potential per particle follows from our theorems about the Vlasov 
limit of the system- specific Boltzmann entropy per particle. 

Next we notice that also the familiar "minimum free energy principle" for 
follows from combining the Legendre-Fenchel transform fll69p with 
our "maximum entropy principle" in Theorem 121 Thus, for the system-specific 
Helmholtz T-potential per particle we find the variational principle 

-^0,(^9) = inf^,(/), (171) 

with 01 = {f e {^Ua n £^ n £^ in i2i)(M3 x A)} the admissible trial densities, 
and 

5.(/) = £(/) + ^^b(/), (172) 

the Helmholtz free energy functional of /, where ^b(/) is Boltzmann's H 
function of /, given in ( l25|l . and £(/) is the energy functional given in (126|1 . 
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It also follows directly from our results that !F-d{f) takes its infimum over the 
set 21, and that any minimizer of f^{f) over 21 is of the form 

/i?(p, q) = Mp)p4q), (173) 

where 

a^p) = (27r^)-texp(-^-il \pf), (174) 
while p^{q) now solves the following fixed point equation on q space, 

/a (-i /a ^a(9> q)Pi?(q)d3^) dg 

with > prescribed. 

We remark that the various possible relationships between the set of maxi- 
mizers of the maximum entropy variational principle and the set of minimizers 
of the minimum free energy variational principle have been discussed in great 
detail in |EHT00t ICETT05"] . Note that this can be (and was) done without 
proving that the maximum entropy variational principle characterizes the limit 
points of Boltzmann's Ergode (|2]) proper. 

We also remark that the existence of the system-specific Helmholtz T- 
potential per particle in the Vlasov limit for the canonical ensemble was shown 
previously by various techniques. Sub-additivity arguments, such as those used 
to prove Theorem [1], are used in |Kie93j . The very strategy which we applied 
to prove Theorems [2] and [3l which not only yields the variational principle for 
the system-specific Boltzmann entropy but also identifies the limit points of 
the sequence of ergodic ensemble measures as convex linear superpositions of 
infinite products of the optimizers for this maximum entropy principle, was 
originally applied in |MeSp82] to the canonical ensemble for Lipschitz continu- 
ous interactions I/^^; subsequently in |Kie93] and in |CLMP92] this approach 
to the canonical ensemble was generalized to less regular interactions including 
the ones studied here; and in |KiSp99] the limit — )■ oo of A^~^ In Z^(jv) (1/-/?) 



was obtained by adapting this strategy (note the different scaling of (3). 
We emphasize that none of these canonical results implies the existence of 
the Vlasov limit for the system-specific Boltzmann entropy per particle, nor 
captures the limit points of the ergodic ensemble measures, unless it is a priori 
known that the ensembles are (convexly) equivalent, i.e. unless it is known 
that e I—)- s^{e) is concave (more on that in section 7). Our results, by contrast, 
hold irrespective of whether e i— )■ s^{e) is concave or not. 

6.3 The Vlasov limit for subergodic ensembles 

Another spin-off, or in this case rather a variation on the theme of our mi- 
crocanonical results is the straightforward generalization of our Theorems to 
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subensembles whose invariant measures are concentrated on sub-manifolds of 
{H = E} determined by further isolating integrals of the Hamiltonian (fT9|) . 
such as angular momentum if the domain A is rotationally symmetric, or 
the Lynden-Bells' invariant |LBLB99l ILBLB04] which occurs in a generaliza- 
tion of the Calogero-Moser model to particles moving in confined by a 
quadratic potential. Hypothesis {H4) does not hold for these interactions, but 
can be replaced by a weaker one at the expense of some extra work. In those 
cases the entropy maximizer factors into a product of a locally (at q) shifted 
Maxwellian on p space and a purely space-dependent Boltzmann factor. The 
shifted Maxwellian which generalizes ( 1271) to include angular momentum is 
known as a "rotating Maxwellian;" in the case of the Lynden-Bells' Hamilto- 
nian one finds a "rotating-dilating Maxwellian." An announcement of these 
results was made in |Kie08j : details will appear in |KiLa09] . 

7 Unfinished business 

In this last section of our paper we point out some open problems related to 
the ones treated here. 

7.1 The maximum interaction entropy principle 

To the best of the author's knowledge, the maximum interaction entropy prin- 
ciple formulated in Proposition |5] is new. As made clear in Theorem |2] it offers 
a way to directly evaluate the usual variational principle of maximum entropy 
with energy constraint. By contrast, the standard approach to evaluate this 
constrained maximum entropy principle has been rather indirect. Namely, 
a Lagrange parameter (basically d) is introduced for the energy constraint, 
yielding the corresponding fix point equation fll75p for the stationary points 
of the free energy functional. After finding all solution families (not just the 
minimizers of the free energy functional), a parameter representation of energy 
and entropy along the various solution families of f ll75p results, among which 
the one with highest entropy for given energy has then to be selected. Clearly 
our new variational approach appears to be more economical than that. 

One of the simplest tasks would be to prove the existence of a unique solu- 
tion to (1251) at sufficiently high energies e. For Coulomb interactions a unique 
solution is expected for all energies, while for (regularized) Newton interac- 
tions multiplicity of solutions is expected for sufficiently low energies. This is 
suggested by the detailed numerical evaluations of the standard principle of 
maximum entropy with constraints for related equations, cf. |SKS951 ICha02j . 
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7.2 Convergence of the ergodic ensemble measures 

We already pointed out in subsection 6.1 that the sequence of ergodic ensem- 
ble measures converges whenever a unique optimizer exists for the maximum 
interaction entropy variational principle in Theorem |2] and Proposition |5l We 
don't see any reason why the sequence of ergodic ensemble measures should 
not converge when the entropy maximizer is not unique, and so we expect 
that the mere existence of limit points concluded in this paper by using weak 
compactness can actually be upgraded to the existence of a hmit. 



7.3 Characterization of the de Finetti— Dynkin measure 

As also noted in subsection 6.1, the decomposition measure u^drlfiir) is a single- 
ton whenever a unique optimizer exists for the maximum interaction entropy 
variational principle in Theorem O In more general situations we have lit- 
tle information on the decomposition measure uldrlfi^), beyond knowing that 
it reduces to ^(dplf)^) and that q{dp\ge) is supported on the maximizers of 
the maximum interaction entropy principle formulated in Proposition [51 Of 
course, we already mentioned earlier that experience with explicitly studied 
physical systems suggests that supp ^(dp|f)e) is either a finite set or a contin- 
uous group orbit of a compact group, but a general proof or disproof seems 
not available. More is known for the canonical ensemble |KuTa84] . and their 
approach should apply to the microcanonical ensemble to determine u^drlfi^). 



7.4 Large deviation principles 

Whenever ^b(/) has a unique minimizer over Sl^, then Theorems [2] and [3] 
imply that 

Prob (^KR [a^U , /f ") > 5) "-^ V(5 > 0, (176) 



where "Prob" refers to the ensemble measure ([2]) with Hamiltonian ( fT9|) . It is 
desirable to improve (11761) to a large deviation principle, a rigorous variation 
on the theme of Einstein's fluctuation formula. Heuristically we expect 

Prob (dKR (a^^I^J^A > 5) X sup e-^(^B(/)-:tfB(/.)) v5 > 0, (177) 

where 2tf = {/ e (qj n ii^ n £Mn £i)(M3 x A) : e(/) = e}\Bs{fe). In 
|EySp93[ lEHTOOl ICETTOS] such a feat was accomplished for the regularized 
microcanonical ensembles at the level of the 1-point functions. The recent ar- 
ticle |EiSch02j establishes some nice large deviation principles for the n-point 
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functions in a strong topology which allows one to handle some singular in- 
teractions. We expect that the conjectured large deviation principle can be 
proved along their lines. 

We also refer to Lanford's article |Lan73] and the books by Varadhan 
|Var84] and Ellis |E1185] for mathematical background on large deviation prin- 
ciples and their applications to statistical mechanics, and to |Tou08] for a more 
recent review. 



7.5 Vlasov limit for the canonical ensemble measures 

Using the very strategy used in this paper to prove our Theorems [2] and [3l 
the Vlasov limit for the canonical ensemble measures associated with f ll57p 
was established in |MeSp82 ICLMP92| IKie93] under various hypotheses on the 



interactions, covering our {H1)-{H5). This raises the question of whether one 
can conclude the convergence of the canonical ensemble measures associated 
with fll57p from the convergence of the microcanonical ensemble measures (or, 
if convergence cannot be shown, the analog for the limit points). Put differ- 
ently, we ask to extend the conclusions reached at the level of the thermody- 
namic functions to the level of the measures. In |EHT00l ICETT05] such a feat 
was accomplished for the canonical ensemble measures in terms of regularized 
microcanonical ensemble measures, using large deviation principle techniques, 
and issues of equivalence of ensembles were addressed. 



7.6 Interactions without lower bound 

By hypothesis (-^2) we allow the pair interactions to diverge when two parti- 
cles approach each other infinitely closely. However, W\{q,q) is only allowed 
to diverge to +oo, which happens with the repulsive Coulomb interactions 
when q ^ q. Divergence of WA{q,q) to —oo is excluded from our analysis, 
because our postulates imply that is bounded below by £g(iV) > —oo. 
In particular, the —oo singularity of the attractive Newton interactions in 
will have to be regularized. 

The canonical ensemble and regularized microcanonical ensembles have 
been controlled under weaker hypotheses, allowing in particular the interac- 
tions to diverge logarithmically to —oo, see |CLMP92| IKie93j for the canonical 
and |CLMP95| IKiLe97[ IKieOOj for the regularized microcanonical ensembles. 
It should be possible to adapt the technical arguments in these papers to estab- 
lish the Vlasov limit for ([2]) for negative logarithmically singular interactions. 



43 



7.7 Unbounded domains 



In |KiSp99| and |ChKiOO] . unbounded A where allowed for the canonical en- 
semble, and our microcanonical theorems should similarly be extendible to 
unbounded domains under a suitable confinement hypothesis which replaces 
hypothesis {H5), presumably 

{H5') Confinement: e"^^'-"^^^ e 2.\Ax A). (178) 

Incidentally, {H5') not only imposes on behavior of U\ as any of its two argu- 
ments is sent to infinity, it also restricts the manner in which [/a can diverge 
to — oo, e.g. when its two arguments approach each other infinitely closely, 
allowing logarithmic divergence. 

7.8 Ergodic ensembles of quasi-particles 

Our analysis does not cover ergodic ensembles of quasi-particle systems hke 
point vortices moving in two dimensions whose Kirchhoff Hamiltonian is of the 
type ([T]) without the sum of |pp terms. The ergodic point vortex ensemble 
measures are of the type 

/if)(d2^X)= (iV!fiV)(£))-'5(£-/f)(X(^)))d2^X, (179) 

A 

where X(^) := (qi,...,q7v) ^ A^, where now A C and d^^X is 2N- 
dimensional Lebesgue measure, and the pair interactions now feature posi- 
tive logarithmic singularities (for a single specie of point vortices). Onsager 
[Ons49] observed that for such systems a critical £ value exists such that the 
map £ H- 5'(£) is decreasing when £ > S,crit, giving rise to negative ensemble 
temperatures. Regularized microcanonical measures for such vortex Hamilto- 
nians have been analyzed in |CLMP95] under an equivalence assumption to 
the canonical ensemble, and in |KiLe97t IKieOO] without such an equivalence 
assumption!^ It is desirable to find a way to handle the proper ergodic en- 
semble for point vortex and other quasi-particle systems for which the sum of 
squares of kinematical momenta is absent from their Hamiltonian, but clearly 
this will require the introduction of new technical ideas. Incidentally, this 
last sentence applies verbatim also to other scalings than Vlasov scaling, in 
particular to the conventional thermodynamic limit scaling explained in the 
introduction. 

25The authors of |CLMP95| use the primitive n^(iv)(£) of ^2'^(«)(£) (i-e. ® with H = I) 
to define a quasi-microcanonical ensemble entropy when £ < tcrit, and for £ > Ecrit they 
use ri,(jv) (oo) — ri,(jv) (£). In [KiLe97, KieOO a Gaussian approximation to S(I — £) is used. 

We also mention |EySp93] where the approximation flj{N} (£) — flj{N) (£ — A£)is used; these 
authors also regularize the logarithmic singularity of the interactions. 
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There is one exception to what we just wrote: precisely at the critical 
energy Ecru of a point vortex system it is a priori known that all the n-point 
measures have densities given by (1/|A|)®". Taking advantage of this fact, 
O'Neil and collaborators |0NR91l IC0N91] found that for a neutral two-species 
system the vicinity of Ecrit oc NhiN can be analyzed directly using 5(J — £); 
it turns out to be a small-entropy regime where S, not S/N, converges to a 
limit when N — > oo, with £ — CA^lnA^ oc A^. Interestingly enough, this scaling 
falls in between the conventional thermodynamic limit and the Vlasov scaling. 
To the author's knowledge, so far these are the only results for point vortices 
obtained for 6{I — £) proper, i.e. without regularization of the Dirac measure. 

Acknowledgment: The author thanks Carlo Lancellotti for his careful read- 
ing of the manuscript and for his comments. This paper was written with 
support from the NSF under grant DMS-0807705. Any opinions expressed in 
this paper are entirely those of the author and not necessarily those of the 
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A Monotonicity of the ground state energy 

In this appendix we will prove two monotonia convergence results about the 
ground state energy which are used in the setup of our construction of the 
Vlasov limit — > oo. The results and their proofs are rather elementary 
and presumably known, and quite likely to be found in the vast literature 
on U statistics; however, my (certainly incomplete) perusal of the pertinent 
literature has not yet met with success@ 
Here is our first proposition. 

Proposition 7. Let A C MP be a bounded and connected domain. Assume 
the following hypotheses regarding U^lq, q): 

(HI) Symmetry: U^iq, q) = U^iq, q) 

{H2) Lower Semi- Continuity: U\{q,q) is l.s.c. on Ax A 

{H3) Sublevel Set Regularity: \®'^[jyUK{q, q) - minf/A < e|) > 

{HA) Local Square Integrability: Up,{q, ■ ) G -2^ {Br{q) n A) V g G A 

where A is normalized Lebesgue measure for A. For N > 2 define the pair- 
specific ground state energy by 



Then the sequence N i— )■ eg{N) so defined is monotonic increasing and con- 
verges to Eg < oo defined by 



Note that Eg as defined in fllSip coincides with Eg as defined in flTB]) when 
D = 3 and is decomposed into the earlier stipulated sum of Va and W\. 



We begin with the mandatory observation that under hypotheses [HI) 
and {H2) the pair-specific ground state energy Eg{N) defined in fllSOp is well- 
defined; i.e. Eg{N) G M (note that {H'i)h{H4) are immaterial here). 

^^In fact, I originally did not expect monotonicity results of the type proved here to hold 
at all. I was prompted to conjecture the results, and then to prove them, by analyzing 
the numerical results of the computations of the (conjectured) ground state energies £g(-/V) 
for Thomson's problem |Tho04| reported in |Aetal97l IPetal97j , which - divided by either 
iV^ or N{N — 1) - arranged themselves monotonically increasing when plotted vs. N . An 
interesting spin-off of the monotonicity of the pair-specific Thomson energies is a necessary 
criterion for minimality which can be used as a test for the empirical numerical experiments. 
After the present paper was submitted I successfully carried out such a test; see [KieQ9b| . 




(180) 




(181) 



Proof of Proposition ^ 
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We next prove the monotonicity of eg{N), with N > 2. Elementary 
(combinatorial) identities and the single inequality that the minimum of a sum 
is not less than the sum of the minima shows that eg{N + 1) > 6g{N), viz. 



SgiN+l) 



mm 

{qiv.Qjv+i} 



1 



{N+1)N 



E 



> 



{N+1)N 



E 



l<fc<Af+l 



l<k<N+l 



mm 



{'Jlv,<7iV+l}\{'Jfe} 



l<i<j<N+l 



N- 



l<i<j<]V+l 



min ^^^f/A(q.,q 



l<i<j<N 



min iv(F3T) EE^A(gi,g,) 



(182) 



and the proof of monotonicity of eg{N) is complete. 

Next, to prove convergence to Eg given by (118 ip we begin by noting that 
under hypotheses (HI), {H2) and (-^4), the ground state energy Eg defined 
in (11811) is well-defined; actually, for this issue we can even relax (HA) to the 
weaker £,1^^{A) condition which is implied by {HA). We now use the density 
of empirical A^-point measures in the weakly compact set of all probability 

2 

measures on A , and the existence (by {H2)Sz{H3)) of a minimizing sequence 
G (2^ (A) for {p,p), to prove convergence Eg{N) 



Eg. We let A^_^)„, denote 



-X. 



(JV) 



the 2-point measure in A' for a ground state Xg^^ = (Oi, q^; Oat, qN)g of A^ 
points in A (which need not be unique), and let A^j^, be any other 2-point 

2 

measure on A with A^ support points. 

We define the linear functional i— )■ U(^p) by 



UCp) = ^|f/A(q,q)^p(d°gd°g). 

Note that for product measures ^p = p*^^ we have 
ll(p^2) = (p, p). 



(183) 



(184) 



Note furthermore that the functional ^p i— )■ 'U(^p) is generally not continuous, 
because we have only (weak) lower semi-continuity of Ua. In particular, while 
any continuous change in the supporting points of the 2-point measure A^(^) 
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2 

on A results in a weakly continuous change of the 2-point measure, the func- 
tional It evaluated at these 2-point measures, i.e. U{A\^[j^^), generally changes 
discontinuously. However, we do have 

e,{N)=U{Af,^,)<U{A%,). (185) 

Now let {pn}neN be a minimizing sequence in (^nC°)(A) for (p, p) = IX(p®^); 
note that it is not necessary to postulate also that Pn — ^ p for any actual 
minimizer p, as this will follow automatically from the proof. Then, by {H3), 
for any e > we can find an such that 'U(p®^) < Sg + e whenever n > n^. 
So pick any e > 0, let n = n^, and let {qk\km be i.i.d. with a-priori measure 
Pn, G (^ n C^6)(A) for each g^. Then by (-^4) the weak law of large numbers 
for U statistics (of order 2) holds |Hoe48] . and so, in probability, 

U(Ag,))'^(Pn.,Pn.><^. + e (186) 

for each e > 0. By (11861) and (11851) we have 

lim sup eg{N) = lim sup U{A^%^ ) < Eg. (187) 

On the other hand, by the compactness of A and the weak* compactness of 
^(A) we can extract a *-weakly convergent subsequence A'"^!jv) ~^ ^ ^(^ )• 

Moreover, since any convergent sequence of n-point measures A''"(iv) necessarily 

^9 

converges to an n-fold product measure, we have ^p = p®^. Now the weak lower 
semi-continuity of XL gives 

liminflI(A(;;,,)>(p,p>>e,. (188) 

Estimates (I187p and (I188p prove convergence Sgi^N) — Eg. 

Convergence and the earlier proved monotonicity of Sg{N) completes 
the proof of Proposition [71 □ 

We notice that our proof of Proposition [7| yields as a "byproduct" that 
(p, p) = Sg. Thus we have the following noteworthy corollary: 

Corollary 2. Any limit point p*^^ of the sequence of ground state 2-point 
measures Nm rninimizes the bilinear form \i{p®'^) = {p,p)- 

Here is our second proposition. 

Proposition 8. Assume the hypotheses on U\{q,q) stated in the previous 
proposition, and in addition assume that Ux> 0- Then the quasi pair-specific 
ground state energy, defined by 

SgiN)^ min ^^^t/^(g„g), (189) 
is a strictly increasing function of N which converges to Eg defined in U81\) . 
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Proof of Proposition O.- 

First of all, eg{N) is as well-defined as eg{N). 

Next, inspection of the monotonicity part of the proof of proposition [7] 
reveals that the same steps as in fll82p now yield 

e^N + l) > (190) 

and eg{N) > because of the here assumed positivity of U\. The strict 
monotonicity oi N ^g{^) ^ow follows because 

{N + 1){N -1) < N^. (191) 

Lastly, since 1 — N^^ — > 1, the limit of eg{N) coincides with that of eg{N). 
This concludes the proof of Proposition [HI □ 



B Decomposition of the finite measures 

Let Qe E ^"(A^) be the weak limit of {g^^l E ^'{A^)}Nm, and let ^{dp\ge) 
be its unique de Finetti-Dynkin-Hewitt-Savage decomposition measure. (If 
{^TV^^jAfGN has several limit points, as accounted for in the main text, the 
following considerations are valid for the associated converging subsequences 
of finite N measures.) We now show that if supp q{dp\gs) is either a finite set or 
a continuous group orbit of a compact group, then for each p E supp ^(dplf)^) 
we can explicitly construct a family of g^^'^lp] E *P*(A^) satisfying 

limV^np]=p^" (192) 
for each n eN, such that for each N eN, 

Q^^i = J f^(^)[p]^(dp|&). (193) 

B.l The support of <^{dp\gi;) is a finite set 

In the simplest case <j(dp|^e) is a singleton, so that gs = pf^, i.e. 

Jim "^5^2^ = pf" V nEN. (194) 

In this case / g''^^[p] ^(dp|f)e) = g^^^lpe] = g^^2^-, and we are done. 

Next, assume that q{dp\g£) is an arithmetic mean of two singletons, viz. 

<^(dp|^e) = ui6p^{dp) + z/25p2(dp) (195) 
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with < z/i = 1 — z/2 < 1, and let dKR^Pi, P2) = D > he the usual 
Kantorovich- Rubinstein distance between pi and p2- Let -Bd/2(pa;) be the KR- 
open ball in *P(A) which is centered at pk and has radius D/2. Now decompose 
= Af u A^, where Ai n A2 = and ^{^^^(Af ) = Uk, such that Af contains 
all points for which ^A^^-* G -B_D/2(Pfc); when is too small there may be no 
such points, but by the weak density in *^^(A) of the empirical one-point mea- 
sures the set of such points G A^ has positive measure when N is large 

enough. In fact, since by hypothesis the weak limit of {q^^J G *P'*(A^)}7vgN is 
given by = u^pf^ + U2pf^ G q3^(A^), it follows that when N 00 then the 
probability w.r.t. g^^J^ that ^A*-^-* G -Bd/2(pa:) approaches z/^. So if we define 



0^''\Pk] = '^,'dxA^ (196) 
and recall that ^)j^|^(A^) = z/^, it follows that 



lim >W[p,] = pf« (197) 

Af— i>oo 

for each n G N and A; = 1 or 2, and such that for each A^ G N, 

0%l = i^iQ^''\pi] + i^20^''\p2l (198) 

which is fll93p in the case that ^ is the arithmetic mean of two singletons. 

The general case of supp <^ being a finite set is treated similarly in an obvious 
manner, with D now the minimum of the set of distances between any pair 
{pk, pi) picked from the support of ^. 



B.2 The support of c,{(ip\Qir) is a continuous group orbit 

For simplicity we assume that we are dealing with a one-parameter continuous 
group G acting on the base space, like S0{2) acting on A; the generalization 
to more comphcated situations (e.g. S'0(3) acting on A) is straightforward. In 
this case we can pick any particular po ^ supp <j and obtain every other (say) 
Pe G supp^ by acting with a group element go E G thusly, pe = Po ° de- The 
de Finetti etc. decomposition of g^ can then be written as an integral w.r.t. 
Haar measure over the group G of the infinite product measures pf^. The 
corresponding finite N presentation is simply obtained by change of variables 
for g^j^2i; through factoring out the group G, which gives each g^'^^ [pg] uniquely. 
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